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1 Introduction 



Since Maldacena's computation of the three-point function produced by an epoch of single- 
field inflation [1], considerable effort has been invested in obtaining correlation functions 
for more complex scenarios. In part this investment has been motivated by the increasing 
sophistication of observational cosmology. But an equally important influence has been the 
drive to understand the implications of our inflationary theories, especially in examples with 
many interacting degrees of freedom. 

The tool used to extract predictions from inflationary scenarios is quantum field theory, 
in which each observable is constructed from correlation functions. Quantum field theory in- 
corporates the effect of fluctuations on all scales, making each correlation function potentially 
sensitive to both short- and long-distance effects. Over very short distances this implies that 
they must be defined carefully by a renormalization prescription, as in flat spacetime. But 
the most interesting outcome of our investment in studying correlation functions has been 
the realization that they have a rich infrared structure which is very different to the case of 
flat spacetime. Some aspects of this structure were reviewed in Refs. [2-4]. 

Nontrivial structure in correlation functions is often associated with the existence of 
kinematic hierarchies. A paradigmatic example from particle physics is a hadronic jet carrying 
some large energy E but small invariant mass mjet- Another is electroweak processes at large 
momentum transfer s ^ M^. (Here s is the Mandelstam variable and Mz is the mass of the 
Z boson.) The kinematic region where ratios such as E/rrijet or s/M^ become large is called 
the Sudakov region. In this region logarithms of large ratios invalidate perturbation theory, 
and to obtain even a qualitative description requires resummation to all orders. 

Analogues of these effects are responsible for the rich infrared structure of inflationary 
correlation functions. Viewed as a sum of Feynman diagrams, a generic n-point function is 
characterized by the (comoving) momenta carried by each propagator. At tree-level these are 
arbitrary linear combinations of the external momenta {ki, . . . , k„}. Another scale is the 
comoving Hubble length r ~ {aH)~^ where a is the scale factor and H = a/a is the Hubble 
rate. Sudakov-like logarithms involving ratios of all these scales appear in each correlation 
function.^ 

When there are no hierarchies among the IqJ and r, each of these scales will be compara- 
ble to a fixed reference scale k^, . If the n-point functions are expressed in terms of background 
quantities evaluated at the horizon-crossing time for fc*, the Sudakov-like logarithms are small 
and the correlation functions are simple. However, if large hierarchies exist it will not be pos- 
sible to find a A;* for which all logarithms are small and the correlation functions develop a 
complex structure. Indeed, we will see that it is necessary to perform an all-orders resumma- 
tion, just as for hadronic jets and high-energy electroweak observables. Understanding the 
details of this enhanced structure is particularly important because correlation functions with 

^Despite the appealing analogy, the logarithms which appear in inflationary correlation functions are not 
Sudakov logarithms in the following sense: in scattering calculations the Kinoshita-Lee-Nauenberg theorem 
implies that sufficiently inclusive observables are infrared safe, due to cancellation between virtual diagrams 
and real soft or collinear emission [5, 6]. Sudakov effects arise in comparatively exclusive observables because 
exclusion of phase space regions prevents complete cancellation. The Sudakov logarithms measure this large 
remainder. In inflation the background flelds themselves have time dependence and there is no meanigful 
sense in which these logarithms arise as a cancellation between real and virtual effects. For this reason we 
shall refer to them as Sudakov-like logarithms. 
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large hierarchies are valuable observables: for example, n-point functions with qi/qj ^ 1 can 
be used to determine nonlinear and stochastic halo bias [7]. But despite their significance, a 
systematic treatment of these effects has not appeared in the literature. 

In this paper we begin a systematic treatment of correlation functions containing Sudakov- 
like logarithms. We focus on logarithms containing the comoving Hubble scale r, which may 
all be rewritten in the form ln|A;=i,r| (perhaps at the cost of introducing other logarithms of 
the form hiqi/k-t). In this form they coincide with the 'time-dependent' logarithms ln|A:=Kr| 
described in Ref. [2]; see also Ref. [4]. We will return to the remaining 'scale' and 'shape' 
logarithms (involving lngj//c* and Ingj/gj, respectively [8]) in a future publication. We work 
to tree-level in quantum effects but to all orders in the slow-roll expansion. In particular, this 
means that masses are accommodated perturbatively. We will comment on the implications 
of this choice as we develop our argument in §§3-5 below. 

Separate-universe formulae. It has been understood for some time that logarithms of 
the form In |A;=Kr| describe time evolution of each n-point function [9-12]. But even though the 
necessity to account for time evolution is well-known, the detailed connexion has remained 
obscure. As a consequence, most discussions of time-dependence are framed in terms of the 
'separate universe' approach to perturbation theory [13-16] — although this has never been 
justified by a derivation beginning from the underlying quantum field theory. One of our 
key conclusions is that any such justification is necessarily a statement about the divergence 
structure of each correlation function to all orders in the slow-roll expansion. 

According to the separate universe prescription, the two- and three-point functions for 
a set of light scalar fields during inflation satisfy [16] 



«(ki)50;3(k2)) 2 T^T^((50i(ki)<5<^,(k2)). (1.1a) 
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(^</.„(ki)5</.^(k2)5</.^(k3)) D |^^|^(,5<A.(ki)<5</.,(k2)<5</.fc(k3)). 



(2vr) 



-|- cyclic permutations. (l-lb) 

We have adopted the notation of Refs. [17, 18], which will be used throughout this paper. On 
the left-hand side, each correlation function is evaluated at some late time t and transforms 
as a tensor in the tangent space associated with the field-space coordinate (j)a{t)- Indices in 
this tangent space are labelled {a, /3, . . .}. On the right-hand side each correlation function 
is evaluated at the horizon-crossing time t^, for the reference scale fc^,, and transforms as a 
tensor in the tangent space associated with the field-space coordinate (/>i(t*). Indices in this 
tangent space are labelled {i, j, . . .}. Where the field-space manifold is curved it is essential to 
preserve this distinction [18]. In this paper we will work with a flat field-space metric except 
in §5.2, but this index convention remains convenient. The bitensors d(f>a/ dcj)* ^ d (pQ^/d(f>^d(p'j 
measure variation of the inflationary trajectory under a change of initial conditions at the 
horizon-crossing time for the reference scale, and the symbol '13' indicates that subleading 
corrections from 'loop' diagrams have been ignored together with disconnected contributions. 

Eqs. (l.la)-(l.lb) make a number of strong assertions. First, Eq. (1.1a) asserts that, 
at all times and accounting for large contributions from any order in the slow-roll expansion. 
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the two-point function can be written as a linear combination of its values at some arbitrary 
earlier time with coefficients d(pa/d(p*. This implies that (at least when decaying modes 
have become negligible) all time-dependent Sudakov-like logarithms can be factorized into 
these coefficients. In fact, we will see that the d(j)a/d(j)* are highly nontrivial time-dependent 
functions and ignoring their contribution would lead to inaccurate conclusions. Here and 
below, 'to all orders in the slow-roll expansion' is a statement about a perturbative expansion 
in terms of slow-roll parameters evaluated at the horizon-crossing time for k-t. 

Second, Eq. (1.1b) asserts that any time-dependent Sudakov-like logarithms appearing 
in the three-point function exhibit a similar factorizable structure. In principle the three-point 
function can have arbitrary dependence on the external momenta ki [1, 19]. Therefore, at each 
order, the Sudakov-like logarithms could enter with different functions of the k^. In contrast 
Eq. (1.1b) asserts that, to all orders in the slow-roll expansion, almost all Sudakov-like loga- 
rithms factorize into the coefficients d(l)a/d(j)* which appear in the two-point function (1.1a). 
Exceptions are allowed only for Sudakov-like logarithms which multiply functions of the ki 
obtainable from a product of two-point functions. For massless fields, or where masses are 
taken into account perturbatively, these are the 'local' momentum combinations (fef^f)"^ 
and its permutations. Even for these exceptions, Eq. (1.1b) requires that the resummation of 
Sudakov-like logarithms produces a result which is related in a nontrivial way to the factor- 
izable coefficients already present in the two-point function. 

Each of these assertions is a complex statement, applicable to all orders in the slow- 
roll expansion, about the structure of infrared divergences which can be produced by the 
underlying quantum field theory. To our knowledge, these statements have not yet been 
justified. The main aim of this paper is to provide such a justification by deducing the time- 
dependent structure of each correlation function directly from the Sudakov-like divergences 
it contains. We will see that this provides a practical and convenient means to compute the 
evolution of each n-point function in cases where the separate universe picture is complicated 
to apply, or an untrustworthy guide for our intuition. 

Summary. In this paper we argue that the Sudakov-like time-dependent logarithms in in- 
flationary correlation functions are qualitatively similar to those occurring in particle physics, 
and therefore can be understood using the same tools. In §2 we sketch the underlying strat- 
egy, focusing on a physical interpretation of the 'factorizable' structure. In §3 we begin to 
develop this argument in detail. The strategy breaks into two parts. 

• First, in §4.2 we prove a factorization theorem which provides constraints on the terms 
in each correlation function which can be logarithmically enhanced. This factorization 
theorem plays the same role as a proof of renormalizability in applying renormalization- 
group arguments to the ultraviolet behaviour of correlation functions: it guarantees that 
all large logarithms can be assembled into a finite number of 'renormalized' functions. 

• Second, in §4.3, we derive a renormalization group equation ('RGe') which can be used 
to determine these unknown 'renormalized' functions. Once the renormalization group 
equations have been solved the full correlation function can be reconstructed. It will 
turn out that the RGBs are equivalent to the transport equations used to determine 
d(f)a/d(j)* and their higher derivatives [17, 18, 20-24]. 
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In §5 we illustrate our method using examples. In §5.1 we discuss the 'matching' problem of 
the separate universe picture, which is often regarded as a means to determine the evolution of 
quantum fluctuations which have become classical after passing outside the horizon. Implicit 
in this point of view is a 'matching' between quantum and classical calculations. Several 
authors have noted that the precise choice of matching time is ambigious, which could lead 
to uncertainties in an accurate calculation [25-27]. In the renormalization-group approach 
there is no need to invoke a matching from quantum to classical evolution, because the entire 
calculation takes place within the framework of quantum field theory.^ We show how this 
gives a unique prescription which resolves the matching ambiguity. 

In §5.2 we illustrate the practical utility of the RGE method by using it to determine 
evolution equations for covariant versions of the coefficients dcpa/dcp* and their higher-order 
generalizations in the presence of a nontrivial field-space metric. These coefficients are not 
straightforward to determine using classical intuition, making the renormalization-group ap- 
proach simple and reliable. This example could be adapted to more complex cases, perhaps 
including the effect of loop corrections, where the usual classical separate-universe arguments 
do not function or are prohibitively difficult to apply. 

In §6 we conclude with a short discussion of our main results. We provide details of the 
next-order calculation of the three-point correlation function (which is required to deduce an 
appropriate renormalization group equation for the coefficient functions) in Appendix A. 

Notation and conventions. Throughout this paper, we adopt units in which c = h = 1. 
Our index conventions for scalar fields are described below Eqs. (l.la)-(l.lb). We work with 
the action 



where V = 1^(0) is an arbitrary potential and Latin indices o, b, . . . , are Lorentz indices 
contracted with the spacetime metric gab- Except where explicitly indicated, scalar field 
indices are contracted using the flat field-space metric (5^^. We will generally set the reduced 
Planck mass Mp = (SvrG)-^/^ to unity. 

2 Why is resummation necessary? 

It is not at all clear from inspection of Eqs. (l.la)-(l.lb) that they involve a resummation 
of time-dependent Sudakov-like logarithms. In this section and the next we explain why re- 
summation is necessary and develop a heuristic approach to these equations. In §4 we will 
consider the same problem from a more formal viewpoint, that of the dynamical renormal- 
ization group. 

To quadratic order, the action for fluctuations of the scalar fields (pa in (1-2), measured 
on uniform-expansion hypersurfaces, can be written 



^Note that this is purely a statement about the calculation of correlation functions. The need to understand 
how fluctuations decohere (presumably solving the Schrodinger's cat problem) remains. 




(1.2) 
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where an overdot denotes differentiation witli respect to cosmic time. Tlie mass-mixing matrix 
iTia/s can be written explicitly in terms of slow-roll parameters [15, 28, 29] 

1 d /q3 \ 2 ■ ■■ 

ma/S = ~ ( J^4>a(t>p\ = Vap - 'i(t)a(t>P " £000/3 - J^(l){a4'(5), (2-2) 

where e = —H/H^ is the usual slow-roll parameter. Inflation occurs whenever e < 1. At low- 
est order in slow-roll the mass-matrix can be related to the 'expansion tensor' introduced 
in Ref. [17], 

Two-point function. It is convenient to work in conformal time, defined by r = dt' /a{t'). 
For the two-point function only a single hierarchy can be generated, measured by the ratio of 
k to the comoving Hubble scale aH ~ —t~^. Therefore |A;r| ~ 1 at horizon crossing. Working 
to next-order in the slow-roll expansion and evaluating the two-point function at least a little 
after horizon-crossing we find 

Hi 

'2k3 



{6M^l)SM^2))r = (27r)35(ki +k2); 



X < 



1 + 2e, ( 1 - 7E - In 



2k 

2 - 7E - ln(-A;*r) - In — 



+ 



(2.4) 



where '• • • ' indicates higher-order corrections which we have not computed. We have set 
k = ki = k2to be the common magnitude of the external momenta and dropped contributions 
which decay like positive powers of r. The scale k^, determines a horizon-crossing time around 
which we have chosen to Taylor-expand background time-dependent quantities such as H. 
Evaluation at the horizon-crossing time for A:* is denoted by a subscript At this stage 
its precise assigment is arbitrary and can be chosen to suit our own convenience. The terms 
involving Ink/k^, and ln(— /^^.r) are our first examples of Sudakov-like logarithms. 

This result has appeared in various forms in the literature. It was given in the single-field 
case by Stewart &; Lyth [30], who worked directly in terms of the conserved comoving-gauge 
curvature perturbation TZ. The corresponding result for field fluctuations in the uniform- 
expansion gauge was given by Nakamura & Stewart [31], who accounted for multiple fields 
and a nontrivial field-space metric. Higher-order corrections were given by Gong &: Stewart, 
who developed an algorithmic approach to compute them using Green's functions of the 
Mukhanov-Sasaki equation [32, 33].'^ 

Callan— Symanzik equation. The reference scale A;* is not physical, and therefore the 
correlation functions cannot depend on it. Invariance under a change in fc* is expressed by 
the Callan-Symanzik equation. Taking all background quantities to be functions of the fields 
(pa, this equation can be written 

^ + ^ A) i^^c, (ki) • • • Ha„ (k„)) = 0. (2.5) 



dlnk^: dlnk^: ^^j^/ 

^Eq. (2.4) of this paper agrees with Eq. (43) of Gong & Stewart [33], although the time-dependent terms 
appear different due to the way Gong & Stewart selected their time of evaluation. 

Results restricted to a two-field model were quoted by Byrnes & Wands [34] and Lalak et al. ]35], who 
worked with an explicit adiabatic-isocurvature basis. Their results neglect the time-dependent term ln(— fc,r) 
and therefore disagree with our Eq. (2.4). Formulae compatible with (2.4) were given by Avgoustidis et al. ]36]. 



-6- 



The differential operator in brackets (• • • ) is simply the total derivative d/dlnk^,. Specializ- 
ing (2.5) to the two-point function yields, to lowest-order in a slow-roll expansion, 

As above, a subscript '*' indicates evaluation at the horizon-crossing point for fc*. Since (2.6) 
holds for any /c^, it expresses the condition e = —H/H^ which is satisfied for any change of the 
background fields (pa which themselves satisfy the background equations of motion. This is 
sometimes expressed by saying that the change should be 'on-shell'. Eq. (2.6) guarantees that 
we are free to change k^, provided that we compensate by adjusting all background quantities 
according to their classical equations of motion. The appearance of classical equations is a 
consequence of our restriction to tree-level processes.^ 

Time evolution. We now return to Eq. (2.4) and consider its behaviour for different values 
of k and r. 

An expression such as (2.4) is said to be computed using fixed order perturbation theory, 
because the calculation is carried to a predetermined order in the slow-roll expansion. When 
k is comparable to off, Eq. (2.6) shows that we can choose A:* to be approximately their 
common magnitude provided we evaluate background quantities near the horizon-crossing 
time for k. Then ln(— fc^jr) and hik/k^, will both be negligible, and a fixed-order expression 
such as Eq. (2.4) is a good approximation. It is dominated by its lowest-order term H'^/2k^. 

More than a few e-folds after horizon-crossing the hierarchy k/aH = \kT\ becomes 
exponentially small. No matter how we choose k^ it is no longer possible to make both 
In /c//c* and ln(— fc^jr) negligible. Therefore we must accept the appearance of large, growing 
contributions which cannot be absorbed into background quantities by a choice of evaluation 
time. Growing terms of this kind are potentially hazardous. They are sometimes described 
as divergences, because on its own ln( — /c*r) becomes unboundedly large as r — )■ 0. This is 
simply an artefact of perturbation theory, in the same way that the Taylor expansion of any 
function 

fix) = f{xo) + f'{xo){x-xo) + --- 

appears to diverge when \x — xo\ ^ 1. In reality, when ln(— /c*t) ~ 1 we must sum an 
infinite number of terms from a fixed-order expression such as (1.1a) before we can determine 
even its qualitative behaviour. Therefore, unless some principle or symmetry enforces a precise 
cancellation, fixed-order perturbation theory can never provide an adequate description. This 
obligation to deal with an infinite sequence of terms of the form [e^, ln(— /c.jr)]"' was emphasized 
by Weinberg [37], who went on to speculate that the methods of the renormalization group 
could be used to perform the resummation, in analogy with the case of Sudakov effects in 
QCD. In this paper we show that this is indeed the case. 

We conclude that any expression which is valid to arbitrarily late times, such as (1.1a), 
must resum an infinite number of terms. For example, consider the well-studied model of 



*One might have expected to recover the equation of motion 3H^a ~ —V,a for each field. This does not 
happen because (2.4) does not depend on the fields individually, but only through their aggregate contribution 
to H. Eq. (2.6) can be regarded as the equation of motion for H. In a model with a single field it is equivalent 
to the classical equation of motion d(j)/dN = — V2e. 
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double quadratic inflation [25, 38-40]. The action is 



(2.7) 



where ip and x ^-re two light scalar fields with m^^ < H. We focus on the dimensionless 



two-point function Safl, defined by 



2P ' 



(2. 



but similar remarks apply to higher n-point functions. 

In Fig. 1 we present a typical super-horizon evolution of E<^(^ calculated using the for- 
mula (1.1a). There is a clearly defined peak in the vicinity oi N = — ln(— A;*r) ~ 20, caused 
by a turn of the trajectory in field space [41], after which asymptotes to a constant. 
Inflation is ongoing through the entire evolution, with e growing to a maximum value ~ 0.26 
near the peak. Fig. 1 shows clearly that Sq,^ can exhibit very rapid evolution. The linear 
approximation involving ln(— A:^,r) is already poor for A'^ > 10 and is totally incorrect (even 
qualitatively) for N > 20. If the initial time had been chosen closer to the peak, the linear 
approximation could easily have been invalidated within < 1 e-fold. 



0.20 - 



0.15 



0.10 - 



0.05 




Figure 1: Super-horizon evolution of the two-point function S^^^ in the double-quadratic model (2.7), normalized by 
the Hubble scale at horizon crossing. This evolution is for the scale which left the horizon (at A'^ = 0) fifty e- 
folds before the end of inflation (at N = 50). The initial conditions were (pi^it = S.2Mp and Xinit = 12.9Mp 
with m^/m^ = 9. 



The flat asymptote for > 20 is associated with convergence to an 'adiabatic limit' 
in which the uniform-density gauge curvature perturbation ^ is conserved [42]. To describe 
this asymptotic plateau certainly requires resummation of an infinite number of powers of 
N = — ln(— /c=kt). In any finite combination the highest power of N must eventually domi- 
nate, leading to power-law growth at large A^. It is only in an infinite sum that sufficiently 
many terms remain available at high-order to balance the growing contribution at low orders, 
enabling T,^^ to be constant at late times. 
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3 Towards a resummation prescription 



We have identified growth of the logarithmic 'divergence' in (2.4) with the infrared region 
where k/aH <C 1. It represents the cumulative effect of interactions which operate over very 
many Hubble times. We would like to separate the physics of these interactions from the 
creation of fluctuations in each /c-mode, which takes place over only a few Hubble times. 
The disparity of these timescales implies that, when they are created, the fluctuations do 
not depend on details of the inflationary model in which they are embedded. It is the sys- 
tematic separation of short-distance, model-independent 'hard' physics and long-distance, 
model-dependent 'soft' effects which produces the factorization in Eqs. (l.la)-(l.lb). 

In this section we illustrate how to perform and interpret factorization of the two-point 
function, and indicate how to extend the method to higher n-point functions. 

3.1 Factorization of the two-point function 

We work with the dimensionless two-point function defined in Eq. (2.8). The for- 

mula (2.4) includes terms at lowest-order in the slow-roll expansion and 0(e) corrections. 
We describe these, respectively, as 'lowest-order' (lo) and 'next-lowest order' (or next-order, 
NLO) terms. Working to NLO, we identify 

5 + 2r*^ - 2<^ ln(-A;.r) - IM^^ In ^) . (3.1) 

The constant matrices and Mq,^ are 

'^a/3 = e*'5a/3(l-7E) + </3(2-7E) (3.2a) 
M*^^ = eJa(S + Kp. (3.2b) 

Heuristic treatment. In what follows, we assume that Sq,^ contains at most powers of 
Ink, but not powers of k itself. This is a consequence of approximate scale invariance. Our 
aim is to separate the long- and short- distance contributions to (3.1). It was explained above 
that when k/aH ^ 1 we cannot make a choice of k^, which eliminates all large logarithms. In 
effect, we do not have enough free scales to absorb every large contribution: it is only possible 
to use k-t to eliminate contributions from one of \n.k/k^ and ln(— A;*r). We choose k^ = k 
to eliminate In/c/fe* and introduce a new, arbitrary 'factorization scale' kp. (However, we 
continue to denote quantities evaluated at the horizon-crossing time for k^, with a subscript 
'*'.) Ultimately, we will see that the factorization scale describes the boundary between long- 
and short-distance physics. We write 



k 

In — -I- ln(— fcpr) 



Hi ( 5ai - <i \n{-k^T) ) ( dpj - u*^j ln{-kpT) ) ( 5ij + 2r*j - 2m*,- In ^ 



(3.3) 



where fa/3 = r^p — MQ^ln2. The second line is valid to next-order, which is the accuracy 
to which we are working. Based on the structure of (3.3), we introduce a 'form-factor' F^i, 
which is a function only of r and the factorization scale kp, 



\p = HlTo^iVp, ( 5ij + 2f *,- - 2m*^- In ^ ) . (3.4) 
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large hierarchy k/ku = \kT\ <C 1 



ky slides towards k ] 



^ ^ \ f 



hard soft increasing wavenumber 



fixed reference scale 'factorization' scale horizon 



Figure 2: Scales appearing in factorization of the two-point function. At late times r — 0, generating a large hierarchy 
between the horizon size kjj = —t~^ and the fixed comoving scale k. The floating scale fcp acts as a 
boundary between the 'hard' (short-distance, k/aH > 1) and 'soft' (long-distance, k/aH < 1) parts of the 
calculation. Beginning with ky = ku, all soft effects are included in the fixed-order part of the calculation. 
The 'renormalization group equation' (3.5) allows us to slide fcp towards k, successively absorbing long- 
distance effects in the Sudakov-like form-factor T^i- Eventually, only 'hard' effects dominated by comoving 
momenta ~ A; are included in the fixed-order part of the calculation. 



We will see that this form-factor modifies the lowest-order functional dependence of 
on background quantities, giving the two-point function an enhanced structure which is not 
visible in fixed-order perturbation theory. Modifications of this type are an inevitable conse- 
quence of large contributions which cannot be absorbed into a change of evaluation point for 
the background quantities already present in a fixed-order expression such as (2.4) or (3.1). 

Since Sq,/? must be independent of fcp it is possible to write a Callan-Symanzik equation 
similar to (2.5). For (3.4) to be independent of requires 

dr„^ _ ^ 

dln/tp" 

In Eqs. (3.3)-(3.5) we have made use of the index convention described in §1. The indices 
a, /?, ... label the tangent space at the late time r, whereas indices i, j, ... label the 
tangent space at the horizon-crossing time for the factorization scale /cp. Note that in (3.5) 
we have rewritten the time of evaluation for Uij as this horizon-crossing time. This is formally 
acceptable because the error we incur is 0{£'^) and therefore below the precision to which 
we are working. In §4 below we will see how to give a more mathematically satisfactory 
formulation of this argument. In this section the discussion is not intended to be rigorous, 
but to highlight the physical meaning of factorization. 

Eq. (3.5) is the 'backwards' evolution equation introduced by Yokoyama et al. [20-22] 
for the separate-universe coefficient d(j)a/dct)i. This discussion shows that it can equally be 
regarded as a kind of renormalization-group equation for the 'Sudakov-like' form-factor r^j. 
Its solution requires an initial condition. This can be deduced from (3.3), which gives T^i = 5ai 
when \kpT\ ~ 1. In combination with (3.5) the initial condition enables us to give a physical 
interpretation of the factorization procedure (see Fig. 2). Our starting point is a fixed-order 
expression such as (3.1). 

We begin with \kpT\ ~ 1. With this choice the form-factors Tai = 6ai are trivial and 
contain no information, whereas the fixed-order piece involving In/c/Zcp accounts for contri- 
butions from both long- and short-distance effects. These are effects arising from the regions 
k/aH > 1 and k/aH < 1, respectively. We now use (3.5) to systematically move kp closer 
to k. This adjustment moves all effects generated up to the point where kp/aH < 1 into 
Fcfj. Effects generated when kp/aH > 1 remain within the fixed-order expression. Since (3.5) 
merely expresses that Sq,^ is independent of kp this does not change its numerical value. We 



(3.5) 



terminate this process when kp ~ k, at which point all soft effects generated when k/aH < 1 
are absorbed into the Sudakov-like form-factor. Conversely, the hard factor in (3.3) receives 
contributions only from times when k/aH ~ 1. Because this hard factor is uncontaminated by 
any other scale, ^ fixed-order perturbation theory can be used to obtain a good approximation. 

Infrared safety. In the special case that the Sudakov-like logarithms vanish order-by-order 
there are no soft effects to absorb in Tai, and Tai ~ dai even for kp ^ k. In these circum- 
stances, the hard factor automatically receives contributions only from k/aH ~ 1. In Ref. [18], 
correlation functions with this property were described as 'infrared safe', by analogy with the 
same situation in scattering calculations. Infrared safe quantities decouple from the complex, 
cumulative interactions of soft modes and depend only on whatever physics dominates the 
hard subprocess. Observables built from such correlation functions have the advantage of 
being theoretically 'clean', but the disadvantage that they probe phenomena operating only 
over a narrow window of wavenumbers. 

This discussion makes clear that the form-factor T^i is the differential coefficient d(j)a/d(j)i, 
and therefore that the factorization scale fep corresponds to the initial time '*' used to con- 
struct the separate universe formulae (l.la)-(l.lb). Its arbitrariness corresponds to the arbi- 
trary location of this initial time-slice. The hard, fixed-order contribution corresponds to the 
two-point function (5(/>"50'^)^,. As we explained above, this receives contributions only from 
times when k ~ aH and describes the hard subprocess by which fluctuations are created. 
The creation process is rapid in comparison with the subsequent evolution, so separation of 
scales implies that its form is universal. All details of this discussion have precise parallels in 
the factorization of Sudakov-like form-factors in QCD [43-46]. 

3.2 Factorization of n-point functions with n ^ 3 

With this background it is possible to consider higher n-point functions. We begin with the 
three-point function. 

Removal of external wavefunctions. For reasons that will become clear below, it is 
helpful to discuss the n-point functions separately from their external wavefunction factors. 
In scattering calculations we sometimes discuss 'amputated' diagrams, meaning that the 
entirety of each external propagator is stripped away. In a nonequilibrium or time-dependent 
setting this cannot be done because half of the propagator participates in a nontrivial integral 
over the temporal position of the vertex to which it is attached. In Appendix A.l we show 
that the Feynman propagator takes the form^ [cf. Eq. (A. 22)] 

(T5c/>.(ki,r)5<^,(k2,r')) = (2vr)3<5(ki +k2) x |<7(^' -)^7/3(^, -0 r< r' ^^ ^^ 

[Wa-y{k,T)w^p{k,T } T <T 

where k = |ki| = |k2| is the common magnitude of the momenta on each external line, the 
operator T denotes time ordering and the matrix-valued mode functions corresponding 



^For this statement to be strictly correct we must choose fc* to remove all other sources of large logarithms, 
in the same way that in scattering calculations we must choose the renormalization scale to do the same. 

^Computation of expectation values requires the 'in-in' or 'Schwinger' formalism, in which all field degrees 
of freedom are doubled and therefore there are four possible time-ordered two point functions. Strictly, 
Eq. (3.6) represents the time-ordered -|— |- two point function. For details, see Appendix A or Refs. [47, 48]. 
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Figure 3: Structure of the three-point function as a 'vertex integral' over individual wavefunction factors. Fig. 3a 
depicts the propagator, which is a product of one conjugated wavefunction (dashed line) and one unconju- 
gated wavefunction (solid line). The cross denotes contraction of these factors by matrix multiplication. In 
Fig. 3b, three 'internal' factors participate in an integral over the vertex whereas the three 'external' factors 
do not. Note that the time at the vertex, r), is always earlier than the external time r. In Fig. 3c we show 
the vertex integral with these external factors removed. 



to each half of the propagator are contracted with each other. In Fig. 3a we depict the 
conjugated mode w*^^ as a dashed half-hne, and the unconjugated mode Wap by a sohd half- 
hne. The index contraction is denoted by a cross joining the dashed and sohd halves. With 
these conventions, the three-point function can be depicted as in Fig. 3b. The external, 
conjugated mode functions (dashed lines) are evaluated at some time r which is taken to be 
later than the time ry associated with the internal vertex. Finally in Fig. 3c we show the 
3-point function with these external, conjugated wavefunction factors removed. The rules of 
the 'in-in' formulation of quantum field theory (see Appendix A) show that the full three- 
point function should be obtained by adding Fig. 3b and its complex conjugate. Therefore 
the three-point function has the structure 

(5(/)a(ki)(50/3(k2)5</'T,(k3))^ = 2Re w*^^{ki,T)wl^{k2,T)w*^{k3,T)Ixf^^{T,ki,k2,k3) , (3.7) 

where the vertex integral /a^jj/ corresponds to Fig. 3c and is determined using the methods 
of the in-in formalism. It represents the cumulative amplitude for three-body interactions 
up to time r, and is given by an integral over the internal wavefunctions (which measure the 
probability of three particles interacting at a point), weighted by geometrical factors (which 
measure the volume in which the interaction can take place) and terms representing the 
detailed structure of each three-body interaction. 

External wavefunction divergences. The discussion in §3.1 shows that although each 
wavefunction Wai3{k,T) contains 'divergent' logarithms as r — )• 0, these are absorbed into 
the form-factor Tai in the combination Ta:j{T, kF)wji{k, —k~^). (Here, the wavefunction Wji 
is evaluated at the horizon-crossing time for the factorization scale fcp.) In concrete terms, 
although the 'divergent' logarithms in Wij{k, —kp^) now appear as powers of Inkp, the com- 
bination Taj{T, kp)'Wji{k, —kp^) is independent of kp. Therefore these 'divergent' logarithms 
cancel. Meanwhile, the overall r-dependence is determined by the solution to the renormal- 
ization group equation (3.5) and does not suffer from the appearance of 'divergent' terms. 
The same applies to w*^ provided we arrange for the phase of Wa/s to be constant at late 
times. 

We now generalize this to the three-point function. Consider the combination 

Taiir, kp)rpj{T, kp)r^k{T, A;p)(50^(ki)50^(k2)50'=(k3))„;,-i (3.8) 
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Figure 4: Form-factors available for the four-point function. 

which appears as one component of the separate-universe formula (1.1b). Comparison with 
Fig. 3b and Eq. (3.7) shows that, by construction, we expect each form-factor to absorb 
the 'divergent' Inkp terms from its corresponding external wavefunction. But if the vertex 
integral 1^^;^ also contributes any factors of In kp then these will not be absorbed. Therefore 
these large, 'divergent' contributions remain, leaving some residual fcp-dependence in (3.8). 

Eq. (1.1b) introduces a new form-factor d'^(pa/0(p*d(/}*, which provides a means by which 
these In kp terms can be absorbed. However, as explained in §1, it imposes stringent conditions 
(to all orders in the slow-roll expansion) on the structure of any Infcp terms which appear 
in the vertex integral — these must all be proportional to the momentum combination 
ki^^2^ (or its permutations) which can be generated by a product of two-point functions. 
Our key task in applying the renormalization-group formalism to the three-point function will 
be to prove a 'factorization theorem' which guarantees that Ix^^, does possess this structure. 

Higher n-point functions. Essentially the same discussion can be given for the four- and 
higher n-point functions, for which the separate universe method gives a structure similar 
to (1.1b). Each n-point function contains a term like (3.8) [49, 50], which absorbs Inkp 
contributions from the external wavefunctions. Once this is done, a finite number of form- 
factors are available to absorb Inkp contributions coming from integration over the vertices 
of the diagram, analogous to the vertex integral I\fj,^. The number of possible form-factors 
typically increases with n, although they are not all independent. 
For example, for the four-point function we can write 



, (ki ) (5^^ (k2 ) 5(?!>T, (ka ) (505 (k4 ) ) 7 



2 Re 



W*axikl,T)wl^{k2,T)w*^{k3, T)w}^{k4„ t)Ix^j,u7t{t, ^1,^2, ^3,^4 



(3.9) 



where the 'vertex' integral Ix^ivtx receives contributions at tree-level from graphs with the 
two topologies shown in Fig. 4. As above, factors of In kp from the external wavefunctions 
of both diagrams can be absorbed by contracting with Foj. Comparison with the formulae 
of Refs. [49, 50] shows that form-factors are available to absorb factors of In/cp from Ix^un 
which are proportional to fc^^/c^^/c;^^ or fej~'^A:^^|ki -|- ksl"^ (or their permutations). These 
two possibilities correspond to the contact and exchange graphs, respectively. However, the 
form-factor for fej~^/c^^|ki -|- ksj"^ is built out of the form-factor used to absorb divergences 
in the three-point function. Only the form-factor for the contact interaction of Fig. 4a is 
'new'. A similar pattern repeats for all higher n, with only the irreducible n-point contact 
graph generating 'new' divergences: all others must be absorbed by form-factors which have 
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already appeared in m-point functions with m < n. It is clear from inspection of Fig. 4b 
and its analogues for higher n-point functions that this arrangement is plausible because all 
diagrams except the contact graph are built out of lower-order vertices whose divergences can 
be described in terms of lower-order form-factors. 

4 The dynamical renormalization group 

In this section we revisit the analysis of §§2-3 from a slightly more formal viewpoint — that 
of the 'dynamical renormalization group'. The original purpose of these methods was to de- 
termine dynamical scaling laws for correlation functions near an out-of-equilibrium critical 
point [51, 52]. However, they have found various applications to inflationary correlation func- 
tions [53-58]. The approach developed here is similar to the discussion of time-dependence 
generated from loop corrections given by Burgess et al. [12]. More recently, Collins et al. sug- 
gested that the dynamical renormalization group method could be used to determine time 
evolution generated by integrating out heavy modes [59]. In addition to its direct applica- 
tion to inflationary models, our analysis demonstrates the steps which would be required to 
achieve this for an arbitrary n-point function. 

Role of 'renormalizability'. To apply the renormalization group requires a guarantee 
that all large logarithms can be absorbed into a finite number of 'renormalized' quantities. 
When applied to ultraviolet behaviour this guarantee is provided either by the criterion of 
renormalizability (in a strictly renormalizable theory), or the fact that only a finite number 
of irrelevant operators need to be kept in order to make predictions at fixed accuracy (in an 
effective field theory). 

In the present case there is no notion of 'renormalizability' and we do not have either 
of these guarantees. The property which replaces renormalizability is factorization, in the 
sense of §§2-3. In particular, what is required is a guarantee that all divergences produced by 
'vertex integrals' such as Ix^^, and Ix^^-k are proportional to a finite number of combinations of 
the external momenta. As we saw in §3.2, the separate universe method makes the assertion 
that all divergences generated by Ixfj.u produce the shape /cj'^A;^'^ (or its permutations), and 
all divergences generated by Ixuun produce the shapes k^^k2^k^'^ or A;^^A;^^|ki + ksl"'^ (or 
their permutations). 

A rigorous proof of these properties, and their analogues for higher n-point functions, 
is a 'factorization theorem'. Theorems of this kind provide the missing guarantee that only 
a finite number of form-factors can be used to absorb all large logarithmically-enhanced 
contributions. They are an integral part of the apparatus used to study perturbative QCD 
processes such as jets and deep inelastic scattering [43-46, 60]. 

4.1 Two-point function 

In this section our purpose is to prove a factorization theorem for the vertex integral Ix^u- 
Before doing so, we briefly return to Eqs. (2.4) and (3.1) and repeat the analysis of §3 in a 
way which can be generalized more easily to the three-point function. 

Ignoring terms which decay like positive powers of |/ct|, the coefficient defined 
in (2.8) depends on k and r only through powers of ln(— /c^kt) and In/c/fc*. Therefore it can 
be interpreted as a Taylor series in these logarithms, expanded around the arbitrary horizon- 
crossing time for A;*. The renormalization group is a method to reverse-engineer a function 
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from the first few terms in its Taylor expansion. Applied to S^j^, this allows an all-orders 
reconstruction from information about the lowest terms in the perturbative series. The term 
'dynamical' merely indicates that the ordinary renormalization-group method is being applied 
to a series expansion in time rather than energy. 

Dynamical renormalization group analysis. To proceed we define a 'hard' contribution 
to Sq,^ which excludes any enhancement due to soft effects from In |A;=i,r|, 



= Hi ( + 2rlp - 2M:^ In ^ ) . (4.1) 



2k\ 
J 

We interpret Sj^^'''^* as a function of k and k^,, but not r. The full two-point function can be 
written 



a/3 - ^a/3 



+ ln(-fe*r), (4.2) 



where '~' indicates that this expression is valid up to NLO accuracy. In particular, we have 
replaced H^dai^ with in the factor multiplying ln(— A;*r). Although there is a potential 
mismatch in this exchange, it would appear at next-next-order and is therefore below the 
accuracy to which we are working. In the solution to the renormalization-group equation this 
will translate to an error below leading- logarithmic accuracy.'' 

We now interpret the right-hand side of (4.2) as a Taylor series expansion around the 
arbitrary horizon-crossing time for k^,. According to Taylor's theorem. 



ln(-A:,r) + • • • , (4.4) 

|fc*T|=l 



and so the coefficient of ln(— /c^,r) is the derivative of evaluated at |/c*r| = 1. Therefore 



dSa/3 



dlnr 



= + S^A nL = Kxn\p + S^.n^- (4.5) 

|fc.T|=l 



A priori this is a statement about the value of the derivative only at the fixed time r = —k~^. 
But precisely because k.^ is arbitrary, the function dSci/3/dlnT||fc^T-|=i obtained in this way 
must be the same function of the horizon-crossing time for A;* as dSf^^/dlnr is of r. Therefore 
we can immediately promote it to a differential equation valid for all r, obtained by removing 
'*' from all terms in (4.5) 

- = SaA'"A/3 + ^p\u\a- (4.6) 

In this equation Sq,^ is a function of k and r, but u^ji is a function of r only. In the literature it 
is often rewritten in terms of the inflationary e-folding time = j H dt using —dlnr = dA^. 



The leading-logarithm approximation resums terms of the form [e* ln(— fc,r)]" for all n, but not 
e[e, ln(— fc,r)]" or smaller. For a renormalization group equation of the schematic form 

- ^ = ^E, (4.3) 
omr 

it can be proved by standard methods (see, eg., Ref. [61]) that resummation of the leading logarithms requires 
knowledge of u to 0(e); next-to- leading logarithms requires knowledge of u to O(e^), and so on. Mathemat- 
ically, our methods will justify the separate universe method to leading-logarithm order. Although this is 
likely to be sufficient for observable inflation, subleading logarithms must become important over very long 
timescales. 
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It might appear that (4.6) could have been obtained by differentiation of Eq. (4.2) with 
respect to r. This procedure would produce an equation which is symbolically identical, 
but where the right-hand side is evaluated at |A;=Kr| = 1. To understand why (4.6) is valid 
for arbitrary r it is necessary to construct (4.6) via comparison with Taylor's theorem, as 
described above. This justifies the more heuristic method used to obtain Eq. (3.5). 

An explicit solution of Eq. (4.6) requires an initial condition. According to (4.4) we can 
obtain this initial condition from the constant term in the Taylor expansion, hard* 8,9 
obtain an accurate estimate using Eq. (4.2) we should set k-t = kto remove large contributions 
from powers of In/c/fc*, which yields^*^ 



where '• • • ' indicates higher-order corrections beyond NLO, which we have not computed, and 
a sub- or superscript 'A:' now indicates evaluation at the horizon-crossing time for k. Together, 
Eqs. (4.2), (4.6) and (4.7) provide an interpretation of the renormalization-group analysis. The 
hard subprocess — here, creation of fluctuations near horizon exit — depends only on the scale 
/c, and therefore fixed-order perturbation theory can be used to obtain an accurate estimate 
such as (4.7). The absence of other scales implies we can be confident that the subleading 
terms represented by '• • • ' are smaller than those we have calculated. This subprocess is 
used as an initial condition for the renormalization-group evolution, which evolves it to lower 
energies — here, represented by the decreasing value of k/aH. This evolution is accomplished 
by successively 'dressing' the hard factor with the cumulative effect of soft processes. 

Eq. (4.6) is the dynamical renormalization group equation for the two-point function. 
It has already appeared in the literature as a 'transport' equation for Sq,^, as described by 
Mulryne et al. [17, 23, 24]. If desired, the methods of Ref. [17] can be used to convert (4.6) 
into an equivalent equation for the form-factor T^i which coincides with (3.5). Therefore, the 
approach of this section is exactly equivalent to that of §3, or the standard 6N formula [16]. 

4.2 A factorization theorem for the bispectrum 

Let us return to the factorization properties of /a^i^. We should interpret divergences in this 
vertex integral to mean that three-body interactions continue to arbitrarily late times, rather 
than being localized to the neighbourhood of horizon-crossing. Therefore the divergences 
originate from a kinematic region where all decaying modes become extinct. In this region the 

^Therefore only the lowest two Taylor terms are required to reconstruct Ec^g. However, these receive 
contributions from all orders in the slow-roll expansion. 

^ In any renormalization-group procedure there is some arbitrariness in setting the initial condition (4.7), 
because we are free to treat some or all of the constant term Sj^^"^'^* as an additive constant rather than the 
zero-order term in the Taylor expansion. This makes no difference in an exact calculation, but can influence 
the result at finite orders in perturbation theory. A prescription for this choice is called a factorization scheme. 
For details, see eg. Ref. [46]. We absorb the entirety of E^^"^"^* as the initial condition. 

^°At |fcr| = 1 there are contributions to Ea/3 from decaying power-law corrections which we have not written 
explicitly. (See also §5.1.) These terms make no contribution to (4.7) because they do not form part of the 
Taylor series at late times. 

The contribution from Ma/s could be removed by choosing k, — 2k. For higher n-point functions, the 
analogous contribution would be removed by choosing fc, = kt. To maximize the accuracy of each initial 
condition one should therefore evolve each n-point function from a marginally different initial time, but for 
the three- and four-point functions this effect will be negligible. 





(4.7) 
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wavefunctions Wq,^ lose their k-dependent information and begin to evolve in the same way. It 
is this property which limits the number of momentum configurations which can be enhanced 
by divergences, and ultimately leads to the possibility of factorization. Mathematically, this 
means that we require only asymptotic information about the dependence of each mode 
function on k and r.^^ 

Asymptotic behaviour. We will prove the factorization property in two parts. The first 
step is to show that each mode function has a 'gap' in its asymptotic expansion in the limit 
\kT\ 0, 

^i{k, r) = iJiik, k,) [Mk, k,,T) + 0(A;r)2] . (4.8) 

Here, ^'j(/c,r) is a mode function for the scale k and may depend on other labels which we 
collectively denote i. The prefactor tpi{k,k^:) is time-independent and adjusted so that each 
Feynman propagator constructed from has the correct normalization. If the calculation 
were carried to all orders in the slow-roll expansion then would be independent of the 
reference scale k^:, but we allow ipi and Ai to include an explicit dependence when truncated 
to any finite order. In addition, ^j(A;, /c*, r) should be slowly varying, by which we mean an 
arbitrary polynomial in Inr. Therefore there is a 'gap' in the asymptotic expansion caused 
by the absence of a term linear in kr. If the elementary wavefunctions satisfy (4.8) then so 
do their cosmic-time derivatives. Where required we treat these as wavefunctions of different 
species, distinguished by the label i. 

For a background spacetime sufficiently close to de Sitter, Eq. (4.8) can be proved using 
the methods of Weinberg's theorem [64]. Weinberg constructed asymptotic approximations 
to the mode function of a massless scalar field by making an expansion in powers of k"^. 
The inclusion of a slowly varying mass promotes the constants in Weinberg's computation 
to arbitrary polynomials of logarithms. Because the expansion is in powers of k"^ there is 
automatically a gap in the asymptotic expansion of the growing mode. If the decaying mode 
were to begin at relative order kr this property could be lost. However, Weinberg's analysis 
shows that the decaying mode begins at relative order (kr)^, leading to the estimate (4.8). 

The same conclusion can be reached by studying the behaviour of solutions to the field 
equation in the asymptotic future r — )■ [65]. The field equation is 

dr^ dr 

Assuming quasi-de Sitter expansion this is an equidimensional equation with solutions (f) ~ r^, 
where 

In the massless case we conclude A = 0, 3. As above, corrections due to a slowly varying 
mass occur as logarithms and subleading terms in the gradient expansion appear as powers 
of {k/a)\ 

^^In models containing extra characteristic scales, such as the horizon scale associated with a step or other 
sharp feature in the potential (see eg., Refs. [62, 63]), this may no longer be true. In particular, it may happen 
that interactions continue outside the horizon until some characteristic time, and then switch off. In these 
models it need not be possible to analyse the divergence structure of the vertex integral using only asymptotic 
information, and extra momentum configurations may receive large enhancements. We do not consider such 
scenarios in this paper. 



+ 2aH-^ +m'^a^(j) = 0. (4.9) 



-17- 



4.2.1 Divergence structure of vertex integral 

The second step uses (4.8) to show that the vertex integral produces divergences only in very 
restricted combinations. This point was emphasized by Weinberg [37], but here we give a 
more detailed analysis. It is possible to regard the theorem proved in this section as a sharper 
version of Weinberg's theorem — providing control not only over power-law divergences, but 
also the way in which logarithmic divergences appear. On the other hand, Weinberg's theorem 
applies to all orders in the loop expansion whereas our argument applies only to the bispectrum 
at tree-level. It would be of considerable interest to strengthen this result to constrain the 
possible momentum combinations which can be logarithmically enhanced at loop level. 

A generic contribution to the vertex integration will mix three 'internal' wavefunctions, 
evaluated at the time r] of the vertex, and three 'external' wavefunctions which are evaluated 
at the time of observation r. The labels for the internal and external parts need not agree. 
For example, they may be mixed by an off-diagonal mass matrix as in (3.6) or (A. 21), or 
the internal wavefunctions may be differentiated with respect to time. Also, depending on 
the assignment of + and — type vertices, half of the wavefunctions will be conjugated with 
respect to the other half. These details are irrelevant for the purposes of our discussion. We 
denote the external labels {£1,62,63} and the internal labels {«i,i2)^3}- In conclusion, the 
late-time behaviour of each diagram can be deduced from a sum of terms of the form 



n ikm)'>Pim (^m) ^em (^m, t) • • • j 

men, 2,3) 



£{1,2,3} 

^ drj 



-00 V 



A-2n • 



feie2e3iii2i3iv) Ailih,r]) -\ ^i2(^2,f?)H Ai^{k3,r]) + 



(4.11) 



where /eie2e3iii2i3 is also taken to be slowly varying in the sense described above. For per- 
turbatively massless fields n is a nonnegative integer arising from powers of the scale factor 
a ~ — I/t. We begin with four powers of a from the integration measure ^—g. Further 
powers appear in combination with spatial gradients, which each gradient producing a power 
of 1/a. Inverse gradients such as d~'^ produce positive powers of o. The precise value of n 
depends on the net number of spatial gradients carried by the operator giving rise to (4.11). 

More than two spatial gradients: n ^ 2. There are only two choices which produce 
divergences in the vertex integral: n = and n = 1. We will study these in detail below. 
For n ^ 2, the vertex integral in (4.11) is convergent. [Clearly /„ itself may still be divergent 
because of Inr terms from the external wavefunction factors Ae^{km, t).] 

4.2.2 Two spatial gradients 

Now consider the case of two net spatial gradients. The primitive degree of divergence of the 
integral in (4.11) is dr]/r]^ ~ l/r. Therefore decaying parts of the internal wavefunctions 
contribute only to the convergent part of the integral, and decaying parts of the external 
wavefunctions produce contributions which decay at least as fast as r. 

Vertex integral. To proceed we isolate the integral in (4.11), which measures the cumu- 
lative effect of three-body interactions. For n = 1 these interactions are suppressed by a 
net-positive number of spatial gradients ~ k/a, and we should expect them to switch off in 
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the superhorizon limit k/a — )• 0. Hence, we anticipate that the integral produces no diver- 
gences. The late-time limit is determined by 




-^feie2e3hi2isiv)^iiih,v)A2ik2,r])Ai.^{k3,r]) convergent. 



The integral has dimensions of wavenumber or inverse time. Because / and the Ai^ are 
only slowly varying, the only scale available with the appropriate dimensions is l/r itself. 
Therefore the divergent part of the integral must scale like l/r multiplied by a polynomial 
in logarithms. For the interactions encountered in typical inflationary theories, Weinberg's 
theorem guarantees that these 'fast', power-law divergences cancel [64]. 

Therefore, after cancellation of the 'fast' divergences, two-gradient operators with n = 1 
behave in the same way as operators with n ^ 2. The vertex integral is dominated by three- 
body interactions which occur around the time of horizon exit, with negligible contributions 
from interactions occurring much later. Any divergent terms in Ii arise only from Inr- 
dependent terms appearing in the external wavefunction factors A^^. 

4.2.3 Zero spatial gradients 

The remaining case is n = 0, generated by operators with zero net spatial gradients. In this 
case there is no suppression in the limit k/a — )• 0, and therefore no expectation that these 
three-body interactions should switch off at late times. 

Vertex integral. As above, we first focus on the integral term in (4.11). There is no longer 
anything to prevent a 'slow', pure logarithmic divergence in addition to the 'fast' terms 
involving inverse powers of r. However, we show that these slow divergences are generated 
by /c-dependent terms drawn from a single internal wavefunction. 

First, consider the contribution from the growing modes of each internal wavefunction. 
This can be written 



The argument of §4.2.2 shows that this diverges like multiplied by a polynomial in loga- 
rithms. We will discuss these terms when we consider the external wavefunction factors. On 
the other hand, because of the 'gap' in the asymptotic expansion (4.8), the contribution from 
the 0(A;r/) ^-suppressed terms from two different wavefunctions is of the form 



Therefore any pure logarithmic divergences must come from the 0{kr])^ term associated with 
a single internal wavefunction i^- The integral has dimensions of k^, and hence these pure 
logarithms (which do not involve powers of r) must be proportional to fc^, because no other 
scales are available. This is the distinctive 'local' shape which we have been seeking. 





convergent. 



An alternative way to reach the same conclusion is to argue that dri r?-^ In™ r? ~ r'^ for any power m, 
which can be proved by induction. 
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External wavefunctions. In addition to these 'local' pure logarithms, the vertex integral 
will produce fast power-law divergences proportional to and t~^. The 'gap' in (4.8) 
guarantees that there are no terms proportional to t~^. The terms can be ignored. 
When multiplied by the growing modes from each external wavefunction they generate 
contributions which are also of order r~^, and Weinberg's theorem guarantees that they must 
cancel [37].^^ When multiplied by decaying terms in each external wavefunctions they yield 
contributions which decay at least as fast as r, and therefore become negligible at late times. 

The situation is different for the divergences. Weinberg's theorem likewise guaran- 
tees that, when multiplied by the growing modes from each external wavefunction, their net 
contribution must cancel. But it is also possible for these divergences to promote decaying 
terms from the external wavefunctions. As we now explain, these also produce only 'local' 
combinations of the momenta. The argument is similar to that for the internal wavefunctions. 

First, consider the case where the divergence combines with decaying contributions 
from two different external wavefunctions. The net contribution will be r^^ x 0(r2)xO(r2) = 
0(r). Therefore these contributions decay at late times and become negligible. It follows 
that pure logarithmic effects can be generated only by the 0(A;t)^ term from a single external 
wavefunction Cm, and on dimensional grounds will be proportional to k^. There is no need 
to keep track of any contributions generated in the same way because these are also 
controlled by Weinberg's theorem. 

4.2.4 Inverse spatial gradients 

We can not yet conclude that the vertex integration produces only divergences proportional 
to a pure power kf of one of the external momenta, because in Einstein gravity some modes 
of the metric are not dynamical: instead, they are removed by constraints. The process of 
solving these constraints can produce inverse spatial gradients d"^ [1, 19]. Working in ADM 
variables, the metric can be written 

ds^ = -N^ + hij{dx' + N' dt){dx^ + dt), (4.12) 

where N is the lapse function and A^* is the shift vector. Only the three-metric hij carries 
independent degrees of freedom. The lapse and shift are determined by constraints. (See 
Appendix A. 3.) To obtain the third-order action it is only necessary to solve these constraints 
to first order [1, 66]. We set = di'Qi + Pu, where fSu is divergenceless. It can be ignored 
for the purposes of the three-point function because it does not contribute to the third-order 
action. The solution for 'di is [19] 

A IT , 12 , 

-^92^1 = -2rHa + (2^ + ) Scl^a. (4.13) 

In principle the shift vector Aj = di'di can appear in operators with zero net spatial gradients, 
generating enhanced shapes formed from rational functions of the external momenta. (If 
present, these ratios of the external momenta would appear as a prefactor in Eq. (4.11) which 
we have not written explicitly.) Although this outcome is compatible with the conclusions 

^^For this conclusion, we require the stronger argument of Ref. [37] which applies to operators proportional 
to a"* ~ 1/r'' in conformal time. To be clear, we would like to emphasize that Weinberg's theorem excludes 
only 'fast' power-law divergences in correlation functions. As discussed explicitly in Ref. [37], it does not 
restrict the appearance of 'slow' logarithmic effects. 
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of §§4.2.2-4.2.3 these shapes are not local and could not be reproduced by separate-universe 
type formulae. To demonstrate that only form-factors compatible with the separate universe 
method are required we must show that these enhanced non-local shapes are absent. 

Whether ??i can produce divergences in the vertex integral depends on its asymptotic 
behaviour. This is an important issue beyond the renormalization-group approach we are 
developing in this paper, because the decay rate of the shift vector is important in any 
attempt to justify the separate universe method. Weinberg gave an argument using a broken- 
symmetry approach [67, 68]. More recently, Sugiyama, Komatsu &: Futamase showed that 
the Einstein equations require the shift vector to decay at late times during inflation [69] . Up 
to next-order, it can be verified by explicit calculation that (4.13) gives 

4 J-f 

—d'^^i = 0{krjf. (4.14) 

This gives i?i = 0(1). The general analysis of Sugiyama et al. extends this conclusion to all 
orders in the slow-roll expansion. 

Absence of time-dependent terms generated by the shift. This decay rate is not 
sufficiently rapid to prevent the generation of time-dependent terms from arbitrary operators 
involving the shift vector Ni. However, because the decay is exponentially fast in cosmic time, 
it will prevent divergences in any operator which involves two or more powers of d'^^i/a?. 
Only operators linear in d'^'di/a'^ can generate time dependence. It can be shown that the 
second-order action is entirely independent of "i?, and that the third-order action contains no 
terms linear in d'^'di/o?'. Therefore the lapse function gives rise to no time-dependent terms in 
the two- or three-point functions, and no nonlocal momentum configurations can be enhanced 
by divergences. 

4.2.5 Factorization for the three-point function 

In conclusion, to all orders in the slow-roll expansion, each diagram contributing to the three- 
point function with external wavefunctions carrying labels {61,62,63} will schematically be 
of the form 



(4.15) 

u{T)kf + v{T)k^ + 'w{T)kl + finite 



616263112213 



where u, v and w are slowly- varying functions of r (and perhaps also k or A;*) given by 
arbitrary polynomials of logarithms. The 'finite' term is independent of r. The precise form 
of u, V, w and the finite piece depends on the interaction under discussion. The notation 
[■ ■ ■ ]ei62e3iij2i3 indicates that both the divergent and finite pieces produced by the vertex 
integral can depend on the labels for the internal and external wavefunctions. Eq. (4.15) is 
the formal statement of the factorization theorem. 



^''in principle, terms such as NiN^ = dididi'di/a^ can generate time dependence, even though they are 
quadratic in i9i. However, such contractions are not generated in Einstein gravity combined with the scalar 
field theories we are considering. 

^^As described in §4.2.3, the time-dependent 'local' contributions u, v and w may include terms generated 
by promotion of decaying terms in an external wavefunction 'i'em due to divergences in the vertex integral. 
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As for the two-point function, perturbative calculations using quantum field theory will 
not provide us with the functional form of u, v, or w. Instead, they generate terms in the 
Taylor series expansion for these functions based at the arbitrary horizon-crossing time for 
kif. The full functions must be reverse-engineered using the methods of the renormalization 
group. 

Eq. (4.15) clearly exhibits the Sudakov-like enhancements generated by divergent log- 
arithms. In agreement with the discussion of §3.2, we see that large Sudakov-like effects in 
each external wavefunction can be absorbed into the form factor Tai. This accounts for the 
first line in (1.1b). Also, we now see that (4.15) guarantees the vertex integral will produce 
precisely the structure identified in §3.2 as a prerequisite for the second line in (1.1b). 

To develop a renormalization-group approach to the three-point function we will not 
pursue the factorization of this form-factor directly. Instead we follow the approach of §4.1. 
Specializing to a multiple-field model in which the external labels {61,62,63} correspond to 
the different species of scalar fields, we see that the lowest terms in the Taylor expansion 
around the horizon-crossing time for A;* are 



+ [^IV'a(fcn^)|']A;,(fcl)A^^(A;2)^;,(fc3) "^""''^ 



dN 



kfln{-k^T) (4.16) 



-|- permutations, 



where 'permutations' includes the terms generated from the Taylor expansion of t'(r)/c2 and 
w{T)k^ in (4.15). The Taylor coefficient Ua/3 associated with time-dependence of the external 
wavefunction factors is already known. Therefore to reverse-engineer the full functional form 
of la/S'y we require only an estimate of the derivative du\^^/dN\^:. This can be read off from 
the ln(— A;=fT) terms in a next-order computation of the three-point function. We collect the 
details of this calculation in Appendix A. 

4.3 Renormalization group analysis for the three-point function 

We can now use (4.16) to apply the dynamical renormalization group to the three-point 
function. The procedure is very similar to the two-point function analysis in §4.1, although 
potentially complicated because the three-point function depends on three distinct external 
momenta ki, /c2, ^3- In addition to a hierarchy between the scale of these momenta and the 
horizon scale aH, there may now be hierarchies among the external momenta themselves. 

Shape and scale effects. The three-point function enforces momentum conservation ki -|- 
k2 -|- k3 = for the wavevectors kj carried by its external lines. Therefore the kj can be 
regarded as forming a triangle. 



In (4.15) we have redefined these contributions by dividing out a factor of the corresponding external growing 
mode Ae^. This is harmless overall, leading only to a redefinition of the associated source term. No such 
terms have yet been encountered in practical calculations, so there has been no need to keep track of this 
redefinition. The first term of this type would arise from subleading corrections to the V"'S(j)^ vertex, which 
already contributes at next-order. Therefore the effect is neai-next-order in fixed-order perturbation theory 
and would contribute only to the next-leading-logarithm solution of the renormalization group equation. (See 
footnote 7 on p. 15.) 
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In addition to the time-dependent logarithms In |/c=kt| which we have studied in §§2-4.2, 
it is now possible to generate logarithms involving ratios of the form \nki/kt and In/ci/Zc*, 
where kt = ki + k2 + k^ is the perimeter of the triangle. (Logarithms such as In/cj/A:* can be 
converted to In ki/kt + In A;^ / /c* . ) These effects were discussed by Burrage et al. [8, 70], who 
identified them with a response to changes in the size or shape of the momentum triangle. 
When any of these hierarchies become large, resummation of the corresponding logarithms 
will endow the three-point function with a new, enhanced Sudakov-like structure. 

We can regard kt as an analogue of the single scale k in (3.1). Indeed, the term ln2k/k^: 
appearing in (3.1) can be regarded as Inkt/k^,, with kt = ki + k2 the 'perimeter' of the 
momentum 2-gon. If the momentum ki associated with each side of the triangle is not 
too different from kt then there is only one relevant momentum scale and situation is very 
similar to that of the two-point function. This is the 'equilateral' configuration. The converse 
situation occurs if one momentum is very much smaller than kt- In this case there are multiple 
hierarchies and computation of the three-point function becomes more complex. In this paper 
we focus only on the case where ki ~ kt and there is a single hierarchy. We intend to return 
to the question of multiple hierarchies in a future publication. 

External wavefunctions. We divide the analysis into time-dependent terms arising from 
external wavefunction factors, and those arising from the vertex integral. The time-dependence 
associated with external wavefunction factors is 'unsourced', in the sense that new contribu- 
tions are not continuously generated. (One way to regard the time-dependent terms generated 
by external wavefunctions is as a resummation of the infinite sequence of Feynman diagrams 
generated by arbitrary insertions of the mass operator rriaj^ScpaScpp on each external line.) In 
contrast, time-dependence arising from the vertex integral is actively 'sourced' by ongoing 
three-body interactions, as described in §4.2. 

The factorization theorem (4.16) shows how to deal with these unsourced terms. We 
define a 'bispectrum' Ba/s^y by 

{6M^l)SM^2)S(t>y{k3)) = {27TfS{k, + k2 + k3)^||^. (4.17) 

The factor (4J3j/i;f)~^ is generated from products of the ipi which normalize each wave- 
function in (4.8). The factorization theorem shows that the 'sourced' and 'unsourced' pieces 
contribute additively to Ba/3-y and can therefore be considered separately. Therefore we ignore 
all 'sourced' contributions from the vertex integral, which will be dealt with below. 

Excluding the remaining soft effects generated by powers of In |fc*r| from each external 
wavefunction, we obtain a 'hard' component -B^^"^* analogous to (4.1). Because there is 
no hierarchy between the external momenta, the hard component can be computed using 
fixed-order perturbation theory. Its lowest-order contribution was calculated in Ref. [19]. 
Alternatively it may be obtained by extracting the lowest-order terms from Eqs. (A. 42), 
(A. 45) and (A. 49). This contribution is analogous to the lowest-order term H^dap in (4.7). 
In Appendix A. 3 we extend the results of Ref. [19] by computing the three-point function to 
next-order. The subleading contributions from this calculation which are not enhanced by 
soft logarithms are analogous to the next-order term r*^ in (4.7). However, for this discussion 
in this section we will not need an explicit expression for B^^^*. 
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Comparison with (4.16) shows that the unsourced contribution to the bispectrum, which 
we denote B^p^, has the Taylor series expansion 



t) ohard* 



* nhard* i * rjhard* i * rjhard* 



ln{-hT) + ---. (4.18) 

We now apply the argument of §4.1 to obtain the dynamical renormalization-group equation 

- ^Ynr ^ '"'axBxii-i + U/sxBaX'y + u^xBajiX- (4.19) 

A suitable initial condition can be determined from the fixed-order perturbative expression 

for B^^^* . As for the two-point function, the methods of Ref. [17] can be used to show 

that (4.19) is equivalent to absorption of the time-dependent logarithms in the form-factor 
r . 

J- CM- 

Vertex integral. We now return to the 'sourced' contributions generated by the vertex 
integral. Collecting terms from Appendix A. 3, we find 

S^;,^''^ 5 -Ht<p, H-Kt) kf, (4.20) 

i 

where a sub- or superscript '*' denotes evaluation at the time |A;^,t| = 1, and u^js^ is obtained 
by differentiating the expansion tensor Mq,^, given in Eq. (2.3), 



It is symmetric under exchange of any two indices. The superscript 'div, vertex' is a re- 
minder that Eq. (4.20) includes only divergent terms from the vertex integral. We see that 
time-dependent logarithms multiply only the local shapes A;^, ^^'^ ^fi ^ required by the 
factorization theorem. The symmetries of ti^^^ reduce these to the combination kf. 

To apply a renormalization-group analysis, we introduce a form-factor for each shape 
which can be logarithmically enhanced. In this case these are the local combinations, 

B^I^Jj"""^ = aa\i3^kf + ap\a-ikl + a^\apkl. (4.22) 

The form-factors aa|^^ are symmetric under exchange of (3 and 7, but need have no other 
symmetries. The full bispectrum is obtained by adding these sourced contributions to the 
unsourced terms obtained from (4.19). 

As usual, exclusion of all terms enhanced by soft logarithms yields 'hard' components 
for the aQ,|^^, and in the absence of other large hierarchies they can also be computed using 
fixed-order perturbation theory. Comparison of (4.16), (4.20) and (4.22) shows that the a^i^^ 
have a Taylor series expansion 



,hard* 



* hard* I * hard* , * hard* . * v-ihard* v-ihard* 
^aA«A|^7 +'"/3AOa|A7 +^'yX0-a\l3X +^aXfM^XI3 



ln(-A:*r) + --- . (4.23) 



We have replaced by a suitable product of Ti^^'^* which preserves the index symmetries. 
As in §4.1 this may generate a small mismatch at next- next-order, which translates to an 
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error in the solution of the renormalization group which is below leading-logarithm accuracy. 
The renormalization-group prescription gives 

da^jj^^ 

dlnr ^ ^"-*'"^l/37 + ^/3Aaa|A7 + U-yXdallBX + UaX/i^X/^^fif (4.24) 

This is the transport equation for aQ,|^^ derived in Ref. [17], where it was demonstrated that 
its solution generates the form-factor d'^4'a/d(pldcl)*. Combining the sourced and unsourced 
components, we reproduce the anticipated expression (1.1b). 

There is some arbitrariness in choosing an initial condition for a^^/^^. We are free to 
treat any local terms which appear in the lowest-order correlation function -B^^^* either as 
contributions to the unsourced part Ba/B^y, or an initial condition for ciai/j^.^^ Whichever 
choice we make the result is the same, because the results of Ref. [17] show that this initial 
condition contributes additively to the final bispectrum. 

5 Applications 

In this section we briefly discuss two examples which illustrate the utility of the renormaliza- 
tion group framework. 

5.1 Matching between quantum and classical eras 

In §§3-4 we emphasized the separation between 'hard', model-independent processes by which 
inflationary fluctuations are created, and 'soft', model-dependent processes by which they 
evolve. The hard creation process is associated with loss of a decaying mode, causing in- 
terference effects to cease and the fluctuations to behave classically. Usually, the separate 
universe picture is considered as a framework which can be used to evolve these classicalized 
fluctuations [14-16]. 

Power-law corrections. Implicit in this point of view is a matching between the quantum 
and classical parts of the calculation. Beginning with Polarski &; Starobinsky [25, 71], several 
authors have noticed that this leads to a potential ambiguity. The issue was later studied in 
more detail by Leach &; Liddle [26] and has recently been revisited by Nalson et al. [27]. 

The problem can be stated simply. In Eqs. (2.4) and (3.1), decaying power-law correc- 
tions which scale like positive powers of |A;t| have been neglected. These already occur at 
lowest order in slow roll, where the equal-time two-point function at an arbitrary conformal 
time r can be written 

{5M^i)SM^2))r = (2vr)35(ki + k2)^(l + k\^). (5.1) 

The term /c^r^ is unity at horizon crossing but decays exponentially fast outside the horizon, 
leading to evolution of the typical fluctuation amplitude {(f)'^)^^'^ by a factor of ^2 between 
horizon crossing and a few e-folds later. By itself this is not significant because one should 
not think of a measurable classical fluctuation with fixed amplitude until the decaying mode 
is lost, which corresponds to fc^r^ becoming negligible [13]. However, it does demonstrate 
that (from this point of view) classical reasoning can not be used until at least a few e-folds 
outside the horizon. 

^®For clarity, we note that this is not the same ambiguity discussed in footnote 9 on p. 16, which was 
resolved by a choice of factorization scheme. 
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Figure 5: Impact of different matching prescriptions on the final value of the power spectrum of ^ in the double- 
quadratic model (2.7). The quoted value is for the scale which passed outside the horizon with field values 
= 8.2Mp and x = 12.9Mp, labelled = on the figures. In each figure, the short dashed line corresponds 
to beginning the classical calculation at Af = 2 e-folds after horizon exit; the medium dashed line at = 4 e- 
folds; and the long dashed line at A'^ = 6 e-folds. The evolution generated by the dynamical renormalization 
group is the grey solid line. In Fig. .5a, initial conditions are chosen using = H'^S^p, with i?jv the 



value of the Hubble rate at the matching time. In Fig. 5b the matching is performed using S^^g = H^S^p, 
where Hi^ is the value of the Hubble rate at horizon crossing. In Fig. 5c the full nlo initial condition (2.4) 
is used, setting fc« = fc to minimize the Infc/fe* terms. In each plot the inset panel shows a close-up of the 
last four e-folds of evolution. 



Choice of matching surface. If the separate universe method is regarded in classical 
terms we cannot apply it at the moment of horizon exit. But Eqs. (2.4) and (3.1) also show 
that we cannot wait too long before switching to the classical calculation, otherwise this 
quantum initial condition will be invalidated by the growing logarithm ln(— A;*r). At the 
earliest, we could perhaps consider the fluctuations to be approximately classical when they 
are ~ 2 e-folds outside the horizon, making k ~ e-4 ~ 0.018. In Fig. 5 we show the effect 
of different matching prescriptions in the double-quadratic model (2.7), using the fc-scale, 
initial conditions and model parameters of Fig. 1. We plot the observable quantity P(^, which 
can be obtained from Sq,^ by a gauge transformation [15]. Each numerical calculation makes 
use of the slow-roll approximation. 

In Fig. 5a we perform the matching at = 2, = 4 and A^ = 6 e-folds after horizon 
crossing. Initial conditions for the classical evolution are set using the lowest-order approx- 
imation with Tia/s = H'^6ai3, where Hj^ is the value of the Hubble rate on the matching 
surface. The decaying power-law term /c^r^ is ignored in setting these initial conditions.^'' 

In Fig. 5b the matching is also performed at A^ = 2, A^ = 4 and A^ = 6 e-folds after 
horizon crossing, but with the initial condition = H'^dap where H}^ is the Hubble rate at 
the moment of horizon exit. Finally, in Fig. 5c we use the full NLO expression (2.4) to provide 
an initial condition at each matching surface. In each figure, the solid grey line shows the 
evolution computed using the dynamical renormalization group (see below). 

With the full NLO initial condition (Fig. 5c) the discrepancy between the asymptotic 
value of computed using each prescription, and also the DRGE evolution, is a few percent. 

^^In §4.2 we argued that the decaying mode began at O(fcr)^, and therefore this term must originate 
from the growing mode which dominates the subsequent classical evolution. Hence, its contribution could be 
included in a self-consistent calculation. This possibility was suggested by Polarski & Starobinsky [25] and 
Lalak et al. [35]. 

We have explained in §4.1 that decaying terms of the form k^r'^ do not contribute the drge evolution, and 
for this reason we have not retained it when setting initial conditions. 
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In Fig. 5b the discrepancy between different matcliing prescriptions is of order 3%-7%. The 
discrepancy between the = 6 prescription and the DRGE evolution is as large as 10%. In 
Fig. 5a the discrepancies are very large, between 10% and 40%. 

In Fig. 5c, the expression (3.1) is used to provide an estimate of at the matching 
surface. One might have expected this estimate to be accurate because it is computed to 
next-order, and the slow-roll parameters are individually small. Therefore we might also have 
expected the final difference in to be negligible, because the separate universe method is 
supposed to be independent of the time at which initial conditions are set. The discrepancy 
arises because, with matching performed at > 0, the fixed-order expression is used to evolve 

up to e-folds. With the matching performed at = 0, all-orders resummation is used. 
Hence, the non-negligible difference between each line in Fig. 5c measures the discrepancy 
between the resummed and fixed-order calculations. It would be even more significant if 
these first A^ e-folds coincided with stronger evolution, such as that occuring near the spike 
in Fig. 1. 

DRGE analysis. As observations improve, an ambiguity of order 10% will certainly be 
important. Therefore we would like to understand which of these prescriptions, if any, is 
correct. This question cannot be answered within a classical framework. The resolution is 
that there is no requirement to divide the calculation into 'quantum' and 'classical' eras joined 
by a matching condition. Although this is a useful picture which guides our thinking, the 
calculation can be carried out entirely within the framework of quantum field theory and 
yields a unique answer. Indeed, as we have seen in §§2-4, the relevant division is not into 
'quantum' and 'classical' effects, but 'hard' and 'soft' processes. It simply happens that, 
because interference effects are absent in the soft superhorizon evolution, it can be described 
at tree-level using the classical equations of motion. In the quantum field theory calculation 
there is no arbitrary matching between 'quantum' and 'classical' eras, but only the floating 
factorization scale. Changes in the factorization scale leave all correlation functions invariant. 

The renormalization-group approach we have described allows each correlation function 
to be calculated purely using the methods of quantum field theory and therefore provides 
a unique resolution of this matching ambiguity. The prescription can be extracted from 
Eqs. (4.6) and (4.7). Positive powers of kr harmlessly decay and do not contribute to the 
Taylor series in ln(— A;,,r) at late times. When reversing-engineering this Taylor expansion, 
the renormalization group evolution should begin at the expansion point k^, and the initial 
condition should be the constant term in the Taylor series. As we argued in §4.1, to obtain the 
most accurate answer we should evaluate (4.7) at the horizon-crossing time for k. Therefore 
we should simultaneously begin the renormalization-group evolution at this time. Hence the 
prescription provided by the DRGE is to use the initial condition Sq,^ ~ but to begin 

the 'classical' evolution at A^ = 0. If desired, higher-order corrections can be retained to give 
a more accurate estimate of the initial value. This has been done to obtain the grey lines in 
Fig. 5, although the change due to inclusion of NLO terms in the DRGE initial condition is 
only of order 0.3%. 

The prescription provided by the dynamical renormalization group is incompatible with 
classical intuition, but by itself this is not important because no classical phase of evolution 
is being invoked in this calculation. 
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5.2 Nontrivial field-space metric 

Elliston et al. [18] recently calculated the inflationary bispectrum generated by a a-model 
Lagrangian with nontrivial field-space metric. The action is 



S 



(5.2) 



where V is a potential and Gai3 is an arbitrary symmetric matrix which can be interpreted 
as a metric. 

To preserve manifest covariance with respect to G^^ it is helpful to describe fiuctuations 
using a variable which obeys a covariant transformation law [72, 73]. Gong & Tanaka gave a 
prescription for rewriting inflationary perturbation theory in terms of such a variable, which 
they denoted Q", and computed the action for fleld fluctuations to third-order [73]. 

Two-point function. The equal-time two-point function for Q"' had already been obtained 
to leading order by Sasaki &; Stewart [15]. The analogue of (2.4) was computed by Nakamura 
&: Stewart [31]. There is a correction to the mass matrix and a more complex structure due 
to preservation of manifest covariance,^* 



(Q"(ki)g^(k2)) D (27r)3<5(ki + k2)|j;^n",n'5,- 



1 + 2e J 1 - 7E - In ■ 



2k 



+ 2w*^ 



2 - ln(-A;=,T) - In 



2k 



7E 



(5.3) 



In accordance with our usual convention we have neglected decaying power-law corrections. 
Greek indices a, f3, ... label the tangent space at fleld-space coordinate (P°'{t), whereas 
indices i, j, ... label the tangent space at coordinate (/)"(r*), where r* is the horizon-crossing 
time for a reference scale k^:. The expansion tensor Wjj satisfles 



w 



+ 



1 1 D 



3/^2 m^a?dt 



■imn] 



H H 



(5.4) 



The indices i, j imply that (5.4) is evaluated at time t^,. The same applies for Gij. The 
derivative D/dt is 0"Vq, where Vq, is the covariant derivative compatible with G^p and 
R'«/375 is the corresponding Riemann tensor. Finally, the parallel propagator Il^j expresses 
parallel transport along the inflationary trajectory in fleld-space between (/)"(r*) and <j)°'{T), 



Texp 



dr? n 



A/3' 



67 



(5.5) 



Indices a', /?', ... label the tangent space associated with the integration time rj, and T is 
the time-ordering operator. For all details we refer to Elliston et al. [18]. 



^*In Ref. [31] the parallel propagator matrices Il'^i were omitted. Because the metric is determined by 
parallel transport, II"iIl''jG'-' — G"''. Therefore, provided we carefully keep track of the evaluation point 
for the metric, these matrices are not needed explicitly at lowest-order. 



- 28 - 



Three-point function. The equal-time three-point function for was obtained up to 
0{(j)/Hf in Ref. [18]. It is 



(Q°(ki)Q'^(k2)Q^(k3)) 5 (2^)^5(ki + k2 + k3 



where A^^'' transforms as a three-tensor in the tangent space at time and is defined by 

I Ai / 1,-. 1,2 j,2 



(5.6) 
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(5.7) 
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We have used fct = fei -|- A;2 + and = = k\k2 + A:iA;3 -|- ^2^3. As in (5.4), all 

time-dependent quantities are evaluated at the horizon-crossing time for fc^,. 

Eqs. (5.6)-(5.7) are considerably more complicated than their flat field-space counter- 
parts."*^^ Growing terms involving ln(— A;*r) are already present at lowest nontrivial order in 
4>lii. At 0(^/i7)^ a double-logarithmic term involving ln^(— A^^kt) appears. In flat field-space, 
growing terms appear only at 0{(f)/H)^ and double logarithms do not occur until even higher 
orders. In combination with the propagator terms n"j these make the time dependence 
of (5.6) appear extremely complicated. 

Dynamical renormalization group analysis. Eqs. (5.3), (5.6) and (5.7) present a case 
study in the use of renormalization group techniques to extract a Sudakov-like structure from 
logarithmic divergences. 

The Taylor theorem on a curved manifold is 

D 



(x) = n V (x, x') • • • n/' (x, x') 



+ dA^">- 



6\ + 



(5.8) 



where the omitted terms are 0{5)?). To write (5.8) we have assumed a smooth path con- 
necting X and x', with Yl" a'{x-,x') the parallel propagator along this path and D/dA the 
corresponding parametric derivative. As in §4, the Taylor theorem (5.8) is our principal tool 
in the renormalization-group analysis. The path connecting x and x' will be the inflationary 
trajectory, and the parametric derivative along it will be D/dN. 

The argument now parallels §4 precisely. Comparing Eq. (5.3) with (5.8) allows us to 
conclude 

D 



diV 



-fe,T = l 



(5.9) 



-Kt=1 

^^The first line in (5.7) covariantizes the lowest-order flat-space computation given in Ref. [19]. If next-order 
corrections were to be kept they would covariantize the computation described in Appendix A. 
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The arbitrariness of allows this to be promoted to a differential equation. Replacing 
by suitable copies of S^^^, as in §4, we find 

= w%ST^ + w'^^ST". (5.10) 

A similar analysis can be performed for the three-point function. Eq. (5.7) shows that 
all time-dependent logarithms multiply only local momentum shapes, in agreement with the 
factorization theorem of §4.2. (This is not altered by the presence of a nontrivial metric.) 
Repeating the argument of §4.3, we conclude 

= ^Nj'ax\p^ + W;3^a„|A^ + ^■y^a^\i3\ + Wo^^Ea/^E^-^, (5.11) 

where Wq,/3^ satisfies 

I ( 0^ 

w„/37 = V(aW^^) + - I V(„R;3|A/,|^)-^-^ - 4R„(^^)A— I . (5.12) 

Jacobi equation. Elliston et al. deduced equations equivalent to (5.10) and (5.11)-(5.12) 
from the evolution of a covariant connecting vector between two adjacent inflationary tra- 
jectories [18], and demonstrated that (when solved perturbatively) these transport equations 
reproduce the 'divergences' involving ln(— /c=|,r) and ln^(— A:*r) in (5.3) and (5.6)-(5.7). The 
analysis in this section shows that this agreement is not restricted to the lowest powers of 
ln(— A;=kt) but will extend to the leading-logarithm terms at all orders. 

Comparing the analysis of Elliston et al. [18] with the discussion given in this sec- 
tion highlights the advantages and disadvantages of the renormalization-group technique. 
Eqs. (5.10)-(5.12) involve combinations of the Riemann curvature tensor and its derivative 
which one could not expect to obtain by accident. The analysis of Elliston et al. [18] pro- 
vides an explanation for the precise combinations which appear, and shows how they can be 
understood in geometrical terms. On its own, however, it does not demonstrate agreement 
to higher orders in ln(— fc=|,r). This is important, because as we argued in §2 and was demon- 
strated explicitly in Fig. 5c, contributions from all powers of \n{—k^T) are required to obtain 
an accurate prediction. 

In contrast, by combining the dynamical renormalization-group framework with a fac- 
torization theorem, we derived a correct set of transport equations directly from the linear 
\D.{—k^T) 'divergences'. (As we explained in footnote 7 on p. 15, there is no need to take 
special measures to account for the double-logarithmic term.) Therefore the DRGE procedure 
gives an algorithmic method to obtain a correct evolution equation, valid to leading-logarithm 
order. But conversely it gives no information about the physical interpretation of this evolu- 
tion, in the same way that the successful calculation of a negative /3-function in QCD did not 
directly give an interpretation of asymptotic freedom as a consequence of antiscreening due 
to gluonic fluctuations. 

In this case Elliston et al. were able to guess the correct physical interpretation of 
Eqs. (5.10) and (5.11)-(5.12). But it is easy to imagine that, in more complex circumstances 
(perhaps where the effect of time-dependence from loop corrections is included), it would no 
longer be possible to guess the correct evolution equation using simple physical arguments. In 
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such cases the dynamical renormahzation group provides a straightforward and algorithmic 
approach. In canonical models it requires the extra complexity of a next-order calculation 
to obtain the first contributions to the linear term in the Taylor expansion [that is, the term 
proportional to ln(— A;*t)]. For the nontrivial metric, the relevant new terms involving the 
Riemann curvature appear already at lowest order. 

6 Summary and discussion 

The separate universe picture is a powerful method to compute the superhorizon evolution 
of inflationary correlation functions. Our main result is a justification of this method from 
the underlying quantum field theory. We have focused on the bispectrum, where the simplest 
discussion can be given, but it is clear what would be involved in extending the analysis to 
higher n-point functions. 

The main steps in the argument are simple. We first show that the vertex integral obeys 
a 'factorization' principle: to all orders in the slow-roll expansion, time-dependent logarithms 
appear only in combination with a few fixed functions of the external momenta. The second 
step is to derive renormalization-group equations for the coefficients of these fixed functions. 
These equations turn out to be the transport equations for the two- and three-point functions, 
which were shown in Ref. [17] to be equivalent to the separate-universe Taylor expansion intro- 
duced by Lyth & Rodriguez. Therefore, in combination, the factorization principle and RGE 
approach reproduce the separate universe method. Technically our argument demonstrates 
this in the leading-logarithm approximation. 

This result justifies the usual practice of inflationary calculations, but also provides a 
systematic means to go beyond it. Burgess et al. [12] studied time-dependent effects from loop- 
level corrections to the power spectrum, making use of the same dynamical renormahzation 
group methods employed in §4. More recently, Collins et al. [59] studied effective field theories 
in a nonequilibrium or time-dependent setting, finding nontrivial contributions to the time- 
dependence after integrating out a massive field. To perform calculations analogous to these 
for n-point functions with n ^ 3 during an inflationary phase it would be necessary to proceed 
as in §4 — first, proving a suitable factorization theorem; and second, constructing the relevant 
renormalization-group equation. (For the two-point function, factorization is usually a much 
more straightforward question due to the absence of vertex integrals.) 

The factorization property admits a clear physical interpretation which is related to 
the usual division of separate-universe calculations into a quantum contribution, generated 
near horizon-crossing, followed by classical evolution. The renormalization-group analysis 
gives a sharper version of this division. Fluctuations are created as they inflate beyond the 
Hubble length, on a timescale which is rapid compared to their subsequent evolution and 
which is insensitive to any model-dependent parameters. Factorization divides the physics of 
the fluctuations into two parts: a 'hard' model-independent subprocess associated with the 
creation event, and 'soft', model-dependent interactions which subsequently dress the hard 
subprocess. One way to understand why correlation functions factorize into self-contained 
parts associated with these 'hard' and 'soft' processes is that the disparity in timescales 
prevents interference effects between the two. 

Because we work to tree-level, it follows that the soft physics can be described using 
the classical equations of motion. Conversely, had we chosen to compute beyond tree-level. 
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the applicability of classical equations of motion might be lost but the division into hard 
and soft physics would remain. Therefore, as a conceptual tool, the hard/soft division takes 
the place of the usual quantum/classical separation, allowing correlation functions to be 
calculated without invoking any classical phase at all. Naturally, this does not absolve us of 
the obligation to explain why the correlation functions obtained in this way can subsequently 
be re-interpreted as the correlations of a classical stochastic process — that is, why we see 
fluctuations which behave classically. The renormalization-group method has nothing new to 
say about this difficulty. 

In the course of our analysis we treated all fields as perturbatively massless. (That 
is, mass terms are treated as interactions.) This has two important consequences. First, 
all fc-dependence occurs as integer power laws modified by logarithms. For example, k^'^^ ~ 
k^[l + dlnk + 0{6^)]; it is the infinite sum of logarithms which reproduces a power-law scaling. 
This restriction to integer powers of k was helpful in proving the factorization theorem in §4.2, 
which depended on the existence of a 'gap' in the asymptotic expansion of each wavefunction 
in integer powers of k, Eq. (4.8). In our applications, the coefficient 5 is of order a slow- 
roll parameter, so in a strict mathematical interpretation we are restricted to scenarios in 
which the slow-roll parameters are small throughout the evolution. Second, the appearance 
of time-dependent logarithmic 'divergences' of the form ln(— /c^kt) is itself a consequence of 
the massless approximation. Incorporation of masses would provide a cutoff which prevents 
these terms growing unboundedly. 

Finally, it should be remembered that our analysis does not apply to all theories. Scenar- 
ios which contain other distinguished scales need not satisfy the conditions of the factorization 
theorem proved in §4.2. To obtain this theorem it was assumed that asymptotic information 
about the wavefunctions provided all necessary data about the superhorizon regime. Where 
this assumption fails, other combinations of the external momenta could receive significant 
enhancements. This corresponds to the well-known restriction that, in applying separate- 
universe arguments, all relevant scales should first be outside the horizon. 

Acknowledgments 

We would like to thank Joseph EUiston, Richard Holman and Andrew J. ToUey for comments 
on a draft version of this paper. 

MD is financially supported by Fundagao para a Ciencia e a Tecnologia [grant SFRH/BD/ 
44958/2008]. RHR acknowledges financial support from Fundagao para a Ciencia e a Tecnolo- 
gia [grant SFRH/BD/35984/2007] in the early stages of this work. RHR also acknowledges 
the hospitality of the University of Sussex. DS acknowledges support from the Science and 
Technology Facilities Council [grant number ST/1000976/1] and the Leverhulme Trust. DS 
acknowledges that this material is based upon work supported in part by the National Science 
Foundation under Grant No. 1066293 and the hospitality of the Aspen Center for Physics. 
The research leading to these results has received funding from the European Research Council 
under the European Union's Seventh Framework Programme (FP/2007-2013) / ERC Grant 
Agreement No. [308082]. 

DS would like to thank Misao Sasaki for the suggestion (made long ago at the workshop 
'Non-Gaussianity from Inflation' held at DAMTP, Cambridge, in April 2006) that the methods 



- 32 - 



of Refs. [32, 33] could be used to determine the superhorizon evolution of the three-point 
function. 



A Three-point function at next-to-leading order 

In this Appendix we collect the components of the three-point function up to NLO in the 
slow-roll expansion, and briefly sketch some details of its calculation. For full details of the 
'in-in' or 'Schwinger-Keldysh' method used to obtain inflationary correlation functions we 
refer to the hterature [1, 47, 48, 64, 74]. 

A.l Propagator 

We require the second- and third-order action for a collection of light, canonically-normalized 
scalar fields coupled to gravity. The second-order action has already been given in (2.1). To 
obtain the propagator G, defined in conformal time by 

G"^(x,r;x',r') = (5,^"(x, r)5/(x', r')), (A.l) 

we should invert the quadratic structure appearing in the action. That gives 

where a' = da/dr. Because spatial translations and rotations are unbroken by the de Sitter 
background, G can depend only on the relative displacement x — x'. Diagonalizing this by 
passing to Fourier space, 

G"^(r,r';x-x') = / Gf {r,r'y^<--'\ (A.3) 

and introducing a dimensionless time x = kr = —k/aH, we find 

where y = kr' and a' now denotes da/dx. Eq. (A. 4) can be transformed into a Bessel equation 
by making the substitution G'^^ = u"'^ \/—x/a. After expanding a to next-order (for details, 
see Refs. [8, 75]) we finally arrive at an equation which can be solved perturbatively, 

where Mq,^ = e^:5ai3 + u^js was defined in Eq. (3.2b). It plays the role of an effective mass 
term, including contributions from coupling to the background geometry. 
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Homogeneous solution. The technology to solve Eq. (A. 5) perturbatively in Mq,^ was 
developed by Nakamura & Stewart [31] and Gong &: Stewart [32, 33]. We adopt their methods, 
with some small modifications noted below to keep each wavefunction purely imaginary in 
the late-time limit. The first step is to find a solution v^'^ of the homogeneous equation. 



dx^ X dx 



9/4' 



3M„ 



0. 



To do so we work perturbatively in Mq,^, writing v 
Up to 0(M) we find 



al3 



X^ 

+ ^1 + 
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I dx^ X dx 



1 d 



dx^ ~'~ X dx ~'~ 



9/4' 



9/4' 



0, 
3 



(A.6) 
0(71^^). 

(A.7a) 
(A.7b) 



Eq. (A. 7a) has solution in terms of Bessel functions of order 3/2. The linear combina- 
tion with correct boundary conditions to match the Minkowski-space wavefunctions in the 



subhorizon limit k/aH 



oo is the Hankel function. 



(A.8) 



where 1J°^I^ is an arbitrary constant matrix. Eq. (A. 7b) can be solved using the causal Green's 
function r(x,z) constructed by Gong & Stewart [32, 33]. This satisfies 



d^r _^ 1 dr _^ 
dx^ X dx 
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9/4 



x^ 



5{x 



and can be written explicitly in terms of the Hankel functions 



r(x,2;) 



ivrz 



X < z 



H. 



(2). 
3/2 '> 



-z)i7^/^(-x) - i/^J^(-z)F^/^(-x) X> Z 



After integrating over F, the solution v^^ can be written 

a/3 



where vi satisfies [32, 33] 



Vl 



vr (— x)^/2 



2e'^ 



'1 -I- ix)e 



dz 



2iz 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



In writing (A.ll) we have adjusted the position of the indices on M to match on both sides of 
the equation. However, since we are generally working with a flat field-space metric, except 
in §5.2, there is no distinction between co- and contravariant indices. 

The scalar part of the zero-order wavefunction ^[^^^^{—x) has the property that it is 
purely imaginary at late times. 



IT ( — x)^/^ 



+ 0(x). 



(A.13) 
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However, in this limit vi is complex. Assembling (A. 8), (A. 11) and (A. 12) gives 



(-x)3/2 



7E - ln(-2x) - i|l ) . (A.14) 



It is the term involving i7r/2 which is problematic. This arises from the asymptotic expansion 
of the exponential integral appearing in Eq. (A. 12). We will see below that to clearly exhibit 
the cancellation of unwanted divergent terms in the calculation of the three-point functions 
it is helpful to adjust the phase of so that this term is removed. Therefore we redefine 
by a pure phase 

^ ^ exp(ir^)i;^^, (A.15) 

where Q'^^ is a constant real matrix. This redefinition does not change the propagator, which 
depends on the combination • v in which this phase cancels out. To remove the imaginary 
term in Eq. (A.14) we should choose 9°^^ = [tt /2)M°'^ . The net effect is to replace vi in 
Eq. (A. 12) with a redefined function (cf. Refs. [8, 66, 76]) 

Feynman propagator. We now return to Eq. (A. 5) and use the homogeneous solution 



V 



a/3 



to construct the Feynman Green's function u'^l^ ?^ In this section it is helpful to switch 
to matrix notation, with (v)q,^ = v"^^ ^ (u)a/3 = u°^^ and so on. In the absence of the mass 
matrix M the field fluctuations are uncorrelated, which implies that we must choose W = 1. 
Together, the Feynman boundary conditions and hermiticity require 

fvt(x') -X- v(2;) x<x' ^^^^^ 
[v"l'(x) • X • v(x') X > x' 

where X is an Hermitian matrix. We will take it to commute with v and v^. This can be 
justified (with our choice 14 = 1) because Eq. (A. 6) implies that both v and are power series 
in the Hermitian matrix M. The normalization condition in (A. 5) likewise guarantees that X 
is a power series in M. Therefore all these matrices commute. Imposing this normalization 
condition, we conclude 



X 



The matrix in brackets (• • • ) is the 'Wronskian matrix' for v and v^^. Abel's theorem shows 
that it can be written 

W(x') = v(x') . - vt(x') . ^ = W(x'o) exp (- £' = W(x^)^, (A.19) 



■^^Strictly, in the context of the Schwinger method this is the time-ordered ++ correlation function. Its 

complex conjugate is the time-ordered correlation function. The H — and — I- correlation functions can 

be extracted from the two alternatives in Eq. (A. 21), bearing in mind that — fields are later than -I- fields for 
the purpose of time ordering, no matter at what physical time they are evaluated. 
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where Xq is arbitrary and W(xg) is independent of M. With our choice 14 = 1, we find 

W(x') = X 1. (A.20) 

7r(— x'j 

Therefore X is Hermitian and commutes with M, in agreement with our earher argument. 
In conclusion, the Feynman propagator can be written 

^ _ ^ (-^)l/2(_^/)l/2 ^ )'vt(a:0.v(a;) x<x' 

4k a(x)a(x') \vt(x)-v(x') x > x' ' 

Uniform expansion of to next-order. We require the lowest-order and next-order 
contributions in slow-roll to G^. This can be obtained from the solutions Vq^ and v'^^ together 
with next-order corrections to the scale factor a, for which purpose we must Taylor expand 
around an arbitrary reference scale as described in §2. Denoting evaluation at the horizon- 
crossing time for /c* by a subscript we find that G^ behaves like the product of two 
effective wave-function factors w, 



w^(x') • w(x) X < x' 
w"l'(x) • w(x') X > x' 



(A.22) 



where w = wq -1- wi satisfies 



wo(x) = --^H^l - ix)e''^ • 1 (A.23a) 



In^(-x)) -lj(l-ix)e^^'-l 



k 



2e^^ + 1^(1 - irE)e^^ - (1 + ix)e-^^ — eH |. 



(A.23b) 



[This result was used in the main text in Eq. (3.6).] Eqs. (A. 23a) and (A. 23b) are both exact 
in X but perturbative in powers of slow-roll and mass M. 

A. 2 Asymptotic behaviour of wavefunctions 

In practical calculations we require the late-time expansion of the elementary wavefunctions, 
w, in order to compute the contribution from external lines. In the present case it is only 
necessary to compute these up to O(x^). We find 

w(x) = ^L=/7,ei^{l + A, -uln^(-x)-ix(l + A,)+ixuln^(-x) + 0(x)2| , (A.24) 
where 

A = 5*(l-7E-ln— jl + u,(2-7E-ln— j, (A.25) 
and u is defined in Eq. (2.3). 
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A. 3 Three-point function 

To compute the three-point function we require the action to third-order. This calculation is 
well-documented in the literature; see, for examples, Refs. [1, 19]. Fluctuations in the fields 
will mix with the metric, which can be accommodated using the ADM decomposition (4.12) 



-N^ dt^ + hijidx' + N' dt){dx^ + dt), 



(A.26) 



where the lapse N and shift vector A^* are determined by constraint equations which follow 
from invariance under temporal and spatial reparametrizations. At the level of the background 
we have N = 1 and A^' = 0. When perturbations are included there are order-by-order 
corrections, discussed in §4.2.4, 



N 



1 + y^Qn, 

n=l 

oo 

di-dn + solenoidal component. 



(A.27a) 
(A.27b) 



n=l 



The solenoidal part of Ni is important only when computing the action at fourth-order and 
above [77, 78]. (Spatial indices i, j, ... are raised and lowered using the three-metric hij.) 
Also, although in principle there could be contributions to the third-order action from On and 
'dn with n ^ 3, the structure of the constraints in the action guarantees that only ai and i?! 
contribute [1, 66]. Specializing to the spatially flat gauge hij 



a^dij, we find^^ 
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D - 
- 2 



d'^x dt a-^ 



3"- 



-I- ai 



,2 -"J 
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(A.28) 



a" 



AH 

+ ai ( 2(/.°5</.„ + ^d'^di + aUGH^ 



where cp'^ = (j)a4'a and repeated indices in the lower position are summed using the flat 
Euclidean metric 6ij. Eq. (A.28) is exact. It has been obtained using an expansion in the 
amplitudes \S(p\ but does not invoke an expansion in slow-roll parameters. The solution for 
ai is 

(A.29) 



and the solution for t?i was used in (4.13), 
AH 



-214 



2H 



2(/)q,5(/)q -|- 2a 



(A.30) 



■^^Eq. (A.28) contains a term linear in d^di/a^, in apparent contradiction with the conclusions of §4.2.4. 
This discrepancy is only an appearance, however. Eq. (A.28) has been simplified using the constraints, which 
allow us to trade i9i for specific combinations of the other metric and scalar fluctuations. Therefore any 
appearance of nonlocality in (A.28) is fictitious and will disappear in the final answer. See the discussion on 
p. 42. 
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We wish to compute the three-point function up to terms which are next-to-leading 
('NLO') in a fixed-order slow-roll expansion. For this purpose we need only extract the LO and 
NLO terms in (A. 28). The LO contributions were given in Ref. [19], where the corresponding 
three-point function was also obtained. They are 



1 



d-^x dt 



bctdjd ^ 



, (A.31) 



(In Ref. [19] the action was quoted in a different but equivalent form.) The new contributions 
at NLO are 

1 



S^^P^ J d'xdt a'(g^fs,H' 



2H^ 
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d- 



(A.32) 



didjd-^ 



where the coefficient tensors gap-f and qa satisfy 

1 0a Vp-y 



3 



2H H'^ 4 " 



0/3 



(A.33a) 
(A.33b) 



(We are using e = (/)V2i^^ which is accurate at lowest order in the slow-roll approximation.) 
In particular, ^a/j^ is symmetric under the exchange /3 o 7, but not under exchanges in- 
volving a. The three-point function for the Va^^ term was obtained by Falk, Rangarajan h. 
Srednicki [9], and later recomputed by Zaldarriaga [10]. These calculations enable the three- 
point function for ga^-^ to be deduced, but results for the remaining terms in Eq. (A. 32) 
have not yet been obtained. 

A. 3.1 Next-to-leading interactions 

Terms in the NLO action (A. 32) contribute to the three-point function beginning at next-to- 
leading order in the slow-roll expansion. Therefore, we need only retain the lowest-order 
contribution (A. 23a) to the propagator. Also we require only the lowest order estimate 
a = — (ff*r)~^ for the scale factor, where as usual we have introduced an arbitrary reference 
scale which serves as the basis for a Taylor expansion of all background quantities, as 
explained in §2. 

Notation. The full set of three-point functions is lengthy. To simplify the presentation, it is 
helpful to factor out a number of common elements. In particular, we define each three-point 
function in terms of a quantity B°'^'^] cf. (4.17), but here we extract a further common factor 
of H^. This must be symmetrized to produce the final result, 

+ [5 permutations]. 



The three-point functions obtained by Falk et al. and Zaldarriaga contained a small typographical error, 
which was corrected in Ref. [79]. The correct expression is given in Eq. (A. 36). 
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where the permutations are formed by simultaneously exchanging the index and momentum 
pairs (Q,ki), and (7,k3). We also define 



kt = ki+k2 + ks, 

= kikj = kik2 + kik^ + k2k3, 



(A.35a) 
(A.35b) 



i<j 



and set 7e to be the Euler-Mascheroni constant, 7e ~ 0.577. 

We highlight terms involving ln(— /c*r) in green. For vertex and internal wavefunction 
corrections, these contribute to the divergent part of the three-point function (4.20). For 
external wavefunctions, they contribute to (4.18). 
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1^2 



ln(-A;*r) + In — + 7e 



(A.38) 
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([k2-k3]2-fciA;i) 



(A.39) 



A. 3. 2 Leading-order interactions 

The LO operators which appear in (A.31) contribute to the three-point function at both 
lowest-order and next-order. The next-order terms come from use of (A. 23b) to correct the 
wavefunctions associated with each internal or external line, or from subleading terms in the 
Taylor expansion of background quantities such as H and (p"-^ including factors of H which 
appear in the scale factor. 

Notation. These sources of next-order corrections were identified and described by Chen 
et al. [66]. Subsequently, rules for explicit evaluation of the integrals which appear at inter- 
mediate stages of the computation were given by Burrage et al. [8, 76], to which we refer for 
all technical details. For our present purposes we require only the functions 



Ji{k3) 



^3 



ln(l - t?3) 
^(^3 + ln[l-^3]), 



(A. 40a) 
(A.40b) 
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where "i^a is defined by 



2k3 

kt ■ 



(A.40c) 



In writing Eqs. (A.40a)-(A.40c) we have used the convention that, although each expres- 
sion depends on all fcj, only the asymmetrically occurring momentum is written explicitly. 
As described in Refs. [8, 76], the Jniks) approach finite values in the limit t?3 — )• which 
corresponds to either of the (physical) squeezed configurations fci — )• or ^2 — ^ 0. 

To describe contributions from the late-time limit of (A. 23b), which are generated by 
external wavefunction factors, we use Eq. (A. 25) and the combination 



■'a , fa 



2 H 



(A.41) 



As above, the following expressions should be symmetrized over simultaneous exchange of the 
index and momentum pairs (a,ki), (/3,k2) and (7,k3). 

• 6(pj36(f)"6<j)a operator 
LO plus NLO from external wavefunctions 



+ 



+ 



A*p^{k2) + Al^{k^) - 24^1n(-fc,T 



^7 



'Pi' 



NLO from vertex 
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(A.43) 



NLO from internal wavefunctions 
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54>°'dj5(j)adjd '^Scpis operator 
LO plus NLO from external wavefunctions 
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NLO from vertex 



B. 



NLO from internal wavefunctions 
ki-ks [klkl 



kt 
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(A.46) 
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where we have defined 

u3„ * 



Cq^^A — —klU*p^5aX — ^2<^/37^aA 

7a^7A = -2kjk'le^5i3^6ax + klk2kz5axu*p.y + klkl5o,x{£Ji3^ + 2m^^) + kiklk2,5 j3^u*^y^ 
lap-yx = klk'iikl + kl)eJaxSfi^ 

^'afi'yX ~ 3felA:2e*<5aA'^/37 



(A.48) 



5(j)ji{di5(f)a)'^ operator. 

This operator is simplified by our choice of phase convention (A. 16). With the un- 
modified wavefunction (A. 12), fast power-law divergences are generated by the internal 
wavefunctions which cancel against other power-law divergences arising from their ex- 
ternal counterparts. 

With the redefined phase in (A. 16), all power-law divergences cancel locally, in the sense 
that we do not need to group terms from different parts of the calculation. In practice 
this is a considerable convenience. 



LO plus NLO from external wavefunctions 
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NLO from internal wavefunctions 
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where we have defined 



^a/37A = 3/cte*5Q,A'5/37 + 2kt5i3^U*^x + 4A;i(^aA''i^7 
^a/37A = ^I3^{kl{k2 + k3)e*6aX " k2k3{eJaX + 2<a)} 

+ 5„A{(i^^ + k2k3)ej,3j - ki{k2 + k3){ejp^ + 2u*p^)} 
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(A.52) 



L 



Relation to factorization theorem. It may appear that only Eqs. (A. 36) and (A. 44) 

satisfy the factorization theorem proved in §4.2. In Eqs. (A.37)-(A.38) and (A. 47) there are 
divergences which do not appear to multiply local combinations of the external momenta. 
[As explained in §4.2, there is no requirement for time-dependence associated with external 
wavefunctions, as in Eqs. (A. 42), (A. 45) and (A. 49), to appear in local combinations.] This is 
caused by the appearance in (A. 28) of a term linear in d'^'Q\ja}. As noted in footnote 21 on 
p. 37, this term appears because the constraints have been used to simplify the third order 
action. Therefore the apparent nonlocality in this term is cancelled by a delicate combination 
of other operators. Eqs. (A. 37) and (A. 47) are the lowest-order and next-order parts of this 
operator and combine to give a local result when symmetrized. The same is true for (A. 38), 
in which d'^^i/a? appears as a proxy for the operator a\4)a5(j)a which can only produce local 
divergences. We could have avoided this complex series of cancellations by working with an 
unsimplified action, at the expense of greater algebraic complexity overall. 
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Note that vertex integral arising from the ScpaidiScpf^)'^ operator yields no time- dependent 
logarithms. This is because the interaction is suppressed by spatial gradients, in agreement 
with the analysis of §4.2. 
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